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ABSTRACT 

The evolution of the metalhcity of damped Lyman a systems (DLAs) is inves- 
tigated in order to explore several scenarios for the nature of these systems. The 
observational data on chemical abundances of DLAs are analysed with robust statisti- 
cal methods, and the abundances are corrected for dust depletion. The results of this 
analysis are compared to predictions of several classes of chemical evolution models 
describing a variety of scenarios for DLAs: one-zone dwarf galaxy models, multi-zone 
disk models, and chemodynamical models representing dwarf galaxies. In order to 
settle constraints for star formation time scales and metal production in DLAs, we 
compare the observational data on the [a/Fe] and [N/a] ratios to the predictions from 
the models. In DLAs, these ratios are only partially reproduced by the dwarf galaxy 
one-zone model and by the disk model. On the other hand, the chemodynamical model 
for dwarf galaxies reproduces the properties of nearly all DLAs. The connection be- 
tween the gas flow evolution and the star formation rate is the reason for the ability of 
this model in reproducing the range of abundance ratios seen in DLAs. The compari- 
son of the observed [a/Fe] and [N/a] trends to the predictions of the chemodynamical 
model is used to derive the formation epoch of dwarf galaxies. The chemodynamical 
model predicts that dwarf galaxies make a significant contribution to the observed 
total neutral gas density in DLAs, and that this contribution is more important at 
high redshifts (z ^ 2 — 3) . This is consistent with a scenario in which the DLA popu- 
lation is dominated by dwarf galaxies at high redshifts and by disks at lower redshifts. 
The relation between DLAs and Lyman Break Galaxies (LBGs) is investigated with 
chemodynamical models describing LBGs. Our results calls for a smoother progres- 
sion in the evolutionary history of DLAs and LBGs rather than a sharp dichotomy 
between the two populations. LBGs and DLAs may constitute a sequence of increasing 
star formation rate, with the LBGs being systems with typically short star formation 
time scales (~ 10* yr), and the DLAs having slower star formation. We also arise the 
possibility that we could be missing a whole population of high HI density column 
objects, with metallicities intermediate between those of DLAs and LBGs. Finally, we 
discuss the possibility that relying only on the observations of DLAs could lead to an 
underestimate of the metal content of the high redshift Universe. 

Key words: cosmology: observations - galaxies: evolution - galaxies: formation - 
galaxies: ISM - intergalactic medium - quasars: general 



1 INTRODUCTION 

The study of chemical abundances in the QSO absorption 
line systems (QALs) and their evolution with redshift are 
key factors to be considered in our attempts to track the 
history of galaxy formation in the Universe. These systems 
are also essential in tracing the cosmic chemical evolution 
because they sample a variety of environments and are rep- 



resentative systems of the early Universe, being enriched by 
the first stellar generations (Pettini et al. 1997a, 2000a). 

The Damped Lyman a Systems (DLAs) are the QALs 
with the highest neutral hydrogen column densities (N(HI) 
^ lO'^" cm~^). These systems are considered to be the most 
suitable for chemical abundance studies in the high-redshift 
Universe because they are believed to be the progenitors 
of present day galaxies, for a number of reasons. Owing to 
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their very large N(HI) values the DLAs dominate the neutral 
content of the universe (Wolfe et al. 1995; Rao & Turnshek 
2000). Furthermore, the comoving baryonic mass density in 
gas inferred from the DLAs at a = 2 — 3 nearly coincides 
with the current baryonic mass density in stars today (Wolfe 
et al. 1995; Storrie-Lombardi & Wolfe 2000). 

The high density column of these systems allowed abun- 
dance measurements of a number of elements never detected 
before in QALs, as is the case of zinc, which has an impor- 
tant role in DLA studies (on the other hand, strong lines 
of abundant elements, e.g. oxygen, become saturated, mak- 
ing it difficult to measure their abundances). In addition, 
their large N(HI) values require little ionization corrections 
(Viegas 1995, Vladilo et al. 2001). In the last few years, ac- 
curate abundances of several chemical elements have been 
derived for DLAs at high z (Lu et al. 1996, Prochaska & 
Wolfe 1999, Prochaska et al. 2001, Centurion et al. 2000, 
Molaro et al 2001), allowing detailed studies of chemical evo- 
lution in these systems, even though there are uncertainties 
that must be taken into account, such as depletion of some 
refractory elements onto dust grains and selection effects. 

A widespread indicator of the metallicity in DLAs is 
the Zn abundance ([Zn/H]), which is believed to be very 
little or not affected at all by dust depletion. The Zn abun- 
dances range from -2 to (Vladilo 1998, Lauroesch et al. 
1996, Pettini et al. 1994, 1997b), suggesting that DLAs 
have a range of formation epochs or star formation histo- 
ries. The form of the evolution of metallicity with redshift 
is not yet clearly known. Several studies argue that at « < 2 
the mean metallicity of damped systems is approximately 
1/10 - 1/30 solar with a large scatter and that it decreases 
at z ~ 3, even though the picture hero is less certain (Pettini 
et al. 1994, Pettini et al 1997b, Pettini 1999, Lu et al. 1996, 
Prochaska & Wolfe 1999). By contrast, in a more extended 
work, Prochaska & Wolfe (2002) concluded that there is 
no significant evolution in metallicity of the DLAs over the 
z ~ 1.5 — 3.5 range, although there appears to be a marked 
decline at z > 3.5, but only by a factor of a few. 

Abundance ratios would be a better tool in DLAs chem- 
ical evolution studios. Unlike absolute values, abundance ra- 
tios do not strongly depend on particular model parameters, 
but mainly on the stellar nucleosynthesis and the adopted 
IMF. The [a/Fe] and [N/a] ratios can be used to settle con- 
straints on time scales for star formation and metal produc- 
tion. While the a elements are mainly produced in type II 
Supernovae (SNe II) in short time scales, the Fe-peak el- 
ements are produced in type la Supernovae (SNe la) in a 
longer time scale. A slow (long time scale) star formation 
correspond to nearly solar values of [a /Fc] while supra- 
solar values may indicate shorter time scales for the star 
formation, which is dominated by SNe II. The [N/a] ra- 
tio can also indicate the star formation time scale, because 
the ejection of nitrogen (mainly produced in intermediate 
mass stars) into the interstellar medium (ISM) is delayed 
with respect to that of the a elements. There is a contro- 
versy, though, about the observed pattern of the [a/Fe] ratio 
in DLAs. Some authors (Pettini et al. 2000a, Ellison et al. 
2001a) detect an [a/Fc] enhancement similar to that of the 
metal-poor stars of our Gala^xy, although the enhancement 
levels are typically lower in DLAs. In contrast, other works 
find nearly solar [a/Fe] ratios, as expected from the chemical 



evolution of dwarf galaxies (Centurion et al. 2000, Molaro 
et al. 2001). 

The discrepancy between the values of the [a/Fe] ratio 
determined by different authors may be related to the way 
they consider dust depletion. It has long been realized that 
dust depletion is a major problem concerning the chemi- 
cal evolution of DLAs, because the observed abundances 
probe only the gas phase of the DLA ISM. Pettini et al. 
(1994) demonstrated that the majority of DLAs exhibit an 
enhancement of Zn relative to Cr, which they interpreted as 
evidence of dust depletion. In addition, DLA surveys gener- 
ally use QSOs selected from magnitude limited samples, in 
which QSOs that lie behind dusty, and presumably metal- 
rich DLAs would be under- represented (Pei & Fall 1995). 
One must know the depletion fraction of each element in 
each system in order to determine its real abundance pat- 
tern; the observed abundance refers only to the fraction of 
the clement in gaseous form. Several authors have proposed 
different approaches to this problem, but there is not yet 
consensus about the procedure to correct for dust depletion. 
Pettini et al. (2000a) applied an uniform dust depletion cor- 
rection to refractory elements in systems for which dust /gas 
ratio is less than 50 %. Savaglio et al. (1999) compared the 
DLAs depletion patterns to those of our own Galaxy in order 
to elaborate a correction for each system separately. Vladilo 
(1998, 2002b), on the other hand, suggested more elaborated 
methods for dust correction based on the dust chemical com- 
position and intrinsic abundances of two elements (Zn and 
Fe) with different dust depletion properties. Previous works 
(Vladilo 1998) assumed an intrinsic solar Zn/Fe ratio, but 
later works (Vladilo 2002b) considered a possible nucleosyn- 
thetic enhancement of this ratio. Recovering the intrinsic 
abundaces of the ISM is a very important goal, but apply- 
ing dust corrections is a complex task, and, therefore, to taJse 
into account the uncertainties in the dust correction process, 
it is better to consider a variety of dust correction models, 
guided by some sensible assumptions about the chemistry 
of dust production in the ISM, and considering a range of 
nucleosynthetic contributions to the Zn/Fe ratio. 

Constraining the nature of the DLAs is crucial for un- 
derstanding the passage of a gas rich universe to the present 
day universe, where the directly observable baryonic con- 
tent of galaxies is concentrated in stars. The DLAs axe the 
dominant reservoir of neutral baryons at every observable 
epoch and are identified as the progenitors of modern galax- 
ies, due to the correspondence of the inferred baryonic mass 
density of DLAs at 2; ~ 2 and the mass density in stars today 
(Storrie-Lombardi & Wolfe 2000). Since DLAs are the most 
directly observable baryonic mass systems at high redshift, 
their properties are an important constraint for theories of 
galaxy formation. One of the most exciting aspects of the 
DLAs is the fact that they provide a means for detecting 
galactic systems in the early stages of evolution. If DLAs 
arc associated to one particular class of objects — e.g. disks, 
spheroids, galactic halos — their study allow us to follow the 
evolution of this population. 

Nevertheless, the nature of the DLAs remains an in- 
triguing problem. Results from early DLA surveys (Wolfe et 
al. 1986, Lanzctta et al. 1991) suggested that these systems 
could be the progenitors of the present day spirals. This 
hypothesis is supported by the kinematics of metal absorp- 
tion lines exhibiting leading edge asymmetry, which could 
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be explained by large rotating disks (Prochaska & Wolfe 
1997, 1998). However, Ledoux et al. (1998) found that the 
correlation of the broadening of the absorption lines 
with the asymmetry of the lines, expected in the disk in- 
terpretation, holds only for ISV < 150 km s^^, while the 
kinematics of the high A.V DLAs is consistent with ran- 
dom motions, suggesting the merging of sub-systems, as ex- 
pected in CDM models. In fact, the predictions of CDM 
cosmologies are often incompatible with the large disk inter- 
pretation. In hierarchical CDM models (Kauffmann 1996), 
high redshift galaxies should be relatively compact and un- 
derluminous in stark contrast with the presence of massive 
disks at 2 ~ 3. The semi-analytical models of Mo, Mao & 
White (1998), in which the DLAs arc identified to ccntrifu- 
gally supported disks at the centre of dark matter halos, 
predict that most DLAs should be low surface brightness 
disks, with a median Vrot < 100 km s~^, pointing to low 
surface brightness galaxies (LSBs) rather than large disks 
as the main candidates for DLAs. On the other hand, the 
DLA data can be used to constrain semi-analytical models 
of gala^xy formation, which include a variety of physical pro- 
cesses (Kauffmann & Haelnelt 2000, Sommerville, Primack 
& Fabcr 2000). In addition, the CDM models have shown 
that the rotating disk interpretation for DLAs is not unique 
and that the kinematical data can also be reproduced by in- 
falling sub-galactic clumps in collapsing dark matter halos 
(Haehnelt, Steinmetz & Ranch (1998). All these studies are 
consistent with a multi-population or multi-component sce- 
nario for the nature of DLAs. A multi-component scenario 
is suggested by the kinematics of the lines of high ions (CIV 
and SilV), which indicate that the high ion subsystem is 
kinematically distinct from the low ion subsystem, although 
both subsystems are interrelated, as would be the case of a 
disk at the centre of a dark matter halo, filled with ionized 
gas radially infalling toward the disk (Wolfe & Prochaska 
2000a,b). 

The multi-population scenario, in which DLAs trace dif- 
ferent types of galaxies, and not a single galaxy class, re- 
ceives support from imaging studies of low redshift {zabs ^ 
1) DLAs (Le Brun et al. 1997; Pettini et al. 2000a; Turnshek 
et al. 2001), which reveal a variety of morphologies for DLA 
host candidates: spirals, compact dwarfs, low surface bright- 
ness galaxies. The luminosities of these objects also span a 
wide range: from a ~ 0.011/^ LSB (Bowen et al. 2001) to a 
3L^ spiral (Le Brun et al. 1997). At higher redshift, the data 
on DLA hosts are scarce. The three hosts of DLA at 2 
studied by M0ller et al. (2002) are typically sub-L* objects, 
and one of them has an irregular morphology. 

The disagreement between some chemical evolution 
studies of DLAs is also suggestive of a multi-population sce- 
nario for DLAs. The controversy about the enhancement of 
the a elements in DLAs illustrates this point. Some stud- 
ies (Pettini et al. 1997) found that the observed abundance 
pattern resembles those of the Galactic thin and thick disk, 
giving support to the disk scenario for DLAs. On the other 
hand, other authors (Molaro et al. 2001, Centurion et al. 
2000) claimed that chemical enrichment of DLAs is char- 
acterized by an episodic or quiescent star formation regime 
similar to predictions of dwarf galaxy models (Matteucci et 
al. 1997). 

In order to understand the nature of the DLAs, there 
have been several works comparing observations of DLAs 



with the predictions of chemical evolution models. Mat- 
teucci et al. (1997) studied the a/Fe and N/0 ratios ob- 
served in DLAs by means of chemical evolution models for 
starburst galaxies and the solar neighbourhood, and con- 
cluded that the most promising models to explain the ob- 
served abundances are those successfully describing dwarf 
irregular galaxies, which assume short but intense bursts 
of star formation. Jimenez et al. (1999) coupled a disk for- 
mation model to a chemical evolution model to investigate 
low surface brightness galaxies and high surface brightness 
galaxies, and conclude that LSBs could dominate the high 
z DLA population. Prantzos & Boissier (2000) and Hou, 
Boissicr & Prantzos (2001) used chemical evolution models 
for disk evolution and showed how observational constraints 
on the metal column density of DLAs lead naturally to the 
observed mild redshift evolution of the abundances of sev- 
eral elements. Cen et al. (2002) used cosmological hydro- 
dynamic simulations to compute the mctallicity evolution 
of DLAs. In their models, the observed slow evolution of 
the DLA metallicity is explained by the steady conversion 
of the highest metallicities sites to galaxies, thus depleting 
this class, while all the lower metallicity systems show, in- 
dividually, an increase in metallicity. Calura, Matteucci & 
Vladilo (2002) compared the observed chemical abundances 
to predictions of one-zone or multi-zone models of chemical 
evolution, and found that DLAs can be explained either by 
disks of spirals observed at large galactocentric distances, 
or by LMC-likc irregular galaxies, or by starburst dwarf ir- 
regulars observed at difi^erent times after the last burst of 
star formation. Finally, we should mention the SPH sim- 
ulations of Nagamine, Springel & Hernquist (2003), based 
on a conservative entropy formulation, and which include 
feedback by supernovae and galactic winds. The results of 
these simulations agree well with observations at « = 3 of 
the DLA abundance, their total neutral mass density, and 
column density distribution. 

In the present work we compare the observations of 
DLAs with the predictions of three classes of models: i) one- 
zone models for dwarf galaxies; ii) multi-zone models for 
disk galaxies; iii) chemodynamical models for spheroids. The 
comparison with the results of the chemodynamical model 
represents a novelty of the present work, since, differently 
from the pure chemical evolution models, the chemodynam- 
ical models allow to calculate the gas flow in a consistent way 
instead of being introduced ad hoc. The connection between 
gas flows and star formation provides a natural explanation 
for the variety of chemical properties of DLAs. In addition, 
in the dwarf galaxy scenario for DLAs, the bursts occurs 
naturally, and one does not need additional model param- 
eters characterizing the burst. Finally, the chemodynamical 
iiKjdc'l. within the scenario in which the DLAs are dwarf 
galaxies, provides the properties of the DLA as a function 
of the galaxy impact parameter, which is not possible in 
the one-zone model. In this sense, the present work comple- 
ments previous works, which, within the disk scenario for 
disks, study the properties of DLAs as a function of the 
galactocentric distance (Prantzos & Boissier 2001, Hou et 
al. 2001, Calura et al. 2003). 

In section 2, we perform a statistical robust analysis of 
the available chemical abundance data for DLAs, and cor- 
rect the abundances for dust depletion. Section 3 presents 
the classes of chemical evolution models used here. In sec- 
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tion 4, we apply models of distinct classes to examine several 
scenarios for DLAs. In Section 5, we summarise our results 
and, within the framework of the chemodynamical model 
of DLAs, we discuss the connection between gas flows and 
chemical evolution of galaxies, the epoch of galaxy forma- 
tion, the evolution of the neutral gas content in DLAs, the 
metal content of the high-redshift Universe, and the relar 
tion between DLAs and LBGs. In this paper we adopt an 
= 0.3, SIa = 0.7, Ho = 70 km s~^ Mpc~^ cosmology 
(the corresponding age of the universe is 13.47 Gyr). 



2 ANALYSIS OF CHEMICAL ABUNDANCE 
DATA OF DLAS 

In the past few years, the intense effort devoted to the 
spectroscopy of DLAs resulted in a large amount of chemi- 
cal abundances derived for these systems. A compilation of 
chemical abundance measurements in DLAs based on the 
literature through early 2003 is shown in Table 1. All of the 
measurements in the table are from high (echelle) resolution 
obtained with 8-10 m telescopes. We consider abundance 
measurements of the following elements: N; Cr, Fe and Zn 
for the Fc-peak elements; O, Si and S for a-elements. Mg 
is not used as an a-elemcrit due to the difflculty in its de- 
termination and the uncertainties in its chemical evolution. 
Table 1 exhibits the HI density column N(HI) (cm~^) in 
logarithm form, and the chemical abundances in the con- 
ventional notation [X/H]=log(X/H)o6s-log(X/H)0. In this 
paper, the reference solar abundances (X/H)© are the N 
and O photospheric values of Holweger (2001), and the me- 
teoritic values of Grevesse & Sauval (1998) for the other 
elements. Whenever the original paper listed in Table 1 uses 
a solar abundance set different from that above, the chemi- 
cal abundances have been recomputed according to the so- 
lar abundances adopted here. We restrict the HI column 
densities to log N(HI)> 20.0, in order for ionization correc- 
tions to be unimportant. Our DLA lower log N(HI) cut-off 
is somewhat below the classical (Wolfe et al. 1986) limit of 
log N(HI)= 20.3. The present work is concerned with the 
[a/N] and [N/a] ratios as diagnostic tools of chemical evo- 
lution of the DLAs. Therefore, Table 1 includes only DLAs 
with at least one chemical abundance measurement (actual 
determination or limit) of an a-element. Silicon is the most 
readily observed a-element, and the Si II 1808 A transition 
is the most widely line used in the determination of the Si 
abundance in DLAs, since it is usually unsaturated allowing 
one to derive reliable Si II column densities. We then set 
z = 0.85 as the low redshift cut-off of our sample, because, 
at this redshift, the SiIIA1808 line falls longward (at 3345 
A) of the atmospheric cut-off near 3200 A. 

2.1 Statistical Analysis 

The data set of Table 1 is large enough to allow for the use 
of statistical methods in analysing abundances and abun- 
dance ratios in DLAs. Following the approach described by 
Ryan et al. (1996) we make use of robust statistical tools 
in order to define the best fit to the data, employing the 
techniques described by Cleveland & Kleiner (1975). Mea- 
surements for which only upper or lower limits are available 
are not included in the analysis. We obtain three resume 



lines that summarise the trend of the data and are defined 
as follows: MM (midmean) - the average of all the points 
between and including the quartiles of the distribution over 
a given range, LSMM (lower semi-midmean) - the midmean 
for all the points below the median and USMM (upper semi- 
midmean) - the midmean for all the points above the me- 
dian. If the data are scattered about the midmean according 
to a normal distribution, the semi-midmeans are estimates of 
the true quartiles of the distribution. In order to obtain the 
resume linos, two procedures, LOWESS and MSTAT. are ap- 
plied to the data. LOWESS is an interactive procedure that 
smoothes a data distribution using a locally weighted robust 
regression according to the method described by Clcvcnd & 
Devlin (1988). This procedure reveals the trend of the distri- 
bution, while MSTAT analyses the spread along the correla- 
tion. In the LOWESS procedure, local regions around each 
data point are identified according to a smoothing parame- 
ter / (0 < / < 1), which controls the fraction of the data 
included in this region. Ideally one should try to choose a 
value of the smoothing parameter that is as large as possi- 
ble without distorting patterns in the data. Suggested values 
are in the range 0.1 — 0.2; in this paper, we use / = 0.1. We 
refer the interested reader to the recent paper by Carrotta et 
al. (2002), which computes summary lines similar to those 
used here, and which gives some additional references on the 
statistical methods and points to a code for the reader to 
use. 

During the chemical evolution the elements of the Fc- 
peak group may follow a similar abundance pattern, as ob- 
served in the Milk Way metal-poor stars where chromium 
and zinc have abundance values close to iron. On the other 
hand, there is an overabundance of Zn compared to Cr and 
Fc in DLAs: in our sample [Zn/H]modian = —1.04 while 
[Fe/H]^edian = -1.58 and [Cr/H]^edian = -1.36 (Figure 
1). This difference is generally attributed to the presence 
of dust in these systems (Pettini et al. 1994, 1997a), which 
would also affect the abundance of a elements, specially of 
silicon, which is highly depleted in the local ISM (Savage & 
Sembach 1996). Sulphur, on the contrary, is almost unde- 
pleted and could be used as a dust-free tracer of a elements 
in DLAs. 

However, the [S/Si] ratio in DLAs has a different be- 
haviour from that expected in the local ISM: the median 
value of the systems that have absolute measurements of 
this ratio are close to solar ([S/Si] median = 0.01). The statis- 
tical analysis exhibits the same trend: the median curve is 
on the solar value hne, while the LSMM and USMM show 
little dispersion around [S/Si] = 0. This fact suggests that 
either the Si dust depletion is not significant in DLAs or 
that the S nucleosynthesis is a result of a contribution dis- 
tinct from that in the local ISM. If both S and Si were 
perfect O tracers, the [S/Si] ratio would have intrinsic solar 
values and the presence of dust (suggested in DLAs by the 
[Zn/Fe] and [Zn/Cr] ratios) would imply [S/Si] > 0. This is 
not the case, as one can see from Figure 2. DLAs exhibit 
both suprasolar and undersolar values of [S/Si]. Also stars 
of the Galactic thick disk show no sign of the a elements 
overabundance when S and Zn are used as tracers of a and 
Fc-peak elements: stars with —1.2 < [Fc/H] < —0.4 have an 
average [S/Zn] of 0.03 (-0.3 < [S/Zn] < 0.3) (Prochaska et 
al. 2000). Matteucci & Frangois (1989) have already inter- 
preted the result of Diaz (1980) that the sulphur gradients 
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Table 1. Chemical Abundances in Dainyjcd Lyman o: SvHtcms. 
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20. 


.62 + 0. 


.01 


-1.50 + 0.05 


-1.88 + 0.04 


-1.69 + 0.01 


-1.17 + 0.03 


<-1.07 


-0.73+0.06 


-1.67+0.06 


3,10 


Q0454 + 039 





,860 


20. 


.69 + 0. 


.06 


-1.03 + 0.12 


-0.89 + 0.08 


-1.02 + 0.08 


-0.80 + 0.12 


_ 


_ 


_ 


1 


PKS0458-020 


2 


.040 


21. 


.65 + 0. 


.09 


-1.19 + 0.09 


-1.50 + 0.10 


-1.77 + 0.10 


>-1.17 


_ 


_ 


_ 


8 


PKS052S-2505 


2 


.141 


20. 


.70 + 0. 


.08 


<-l,09 


-1,29 + 0,09 


-1.26 + 0.27 


-1.00 + 0.09 


-1,07 + 0.09 


— 


-2.05+0.11 


2,19 


PKS0528-2505 


2 


.811 


21. 


.20 + 0. 


.10 


-0,78 + 0,12 


-1,24 + 0,16 


-1.25 + 0.15 


-0.76 + 0.11 


-0,84 + 0.10 


>-2.00 


<-1.47 


2,19 


Q0551-366 


1 


.962 


20. 


.50 + 0. 


.08 


-0.15 + 0.09 


-0,92 + 0,10 


-0,95 + 0.09 


-0.44 + 0.10 


-0.32 + 0.14 


- 


- 


14 


HS0741+4741 


3 


.017 


20. 


.48 + 0. 


.10 


- 


- 


-1.93 + 0.10 


-1.69 + 0.10 


-1.68 + 0.10 


>-1.46 


-2,44 + 0,10 


8,10 


Q0836 + 11 


2 


.465 


20. 


.58 + 0. 


.08 


<-1.13 


<-1.35 


-1,40 + 0.10 


-1.15 + 0.11 


<-1.12 


>-1.82 


- 


8 


Q0841+129 


2 


.375 


20. 


.95 + 0. 


.09 


-1.51 + 0.10 


-1.55 + 0.10 


-1.58 + 0.10 


-1.27+0.10 


-1.38 + 0.09 


>-1.97 


-2.26 + 0.10 


8,19 


C^0841 + 129 


2 


.476 


20. 


.78 + 0. 


.10 


-1.40±0.10 


-1.61±0.10 


-1.75 + 0.10 


1.39 + 0.14 


1.39+0.10 


>-1.92 


-2.59+0.10 


8,19 


Q0930 + 2858 


3 


.235 


20. 


.18 + 0. 


.10 


- 


- 


- 


>-1.80 


-1.71 + 0.15 


>-2.38 


-2.29+0.15 


4 


BRI0951-0450 


3 


.857 


20. 


.60 + 0. 


.10 


- 


- 


-2.00 + 0.10 


>-1.53 


- 


- 


- 


8 


BRI0951-0450 


4 


.203 


20. 


.40 + 0. 


.10 


- 


- 


<-2.59 


-2.62 + 0.10 


- 


- 


- 


8 


PSS0957+33 


4 


.178 


20. 


.50 + 0. 


.10 


- 


- 


-1.87 + 0.11 


-1.50 + 0.10 


-1.31 + 0.12 


>-1.90 


- 


8 


PSS0957+33 


3 


.279 


20. 


.32 + 0. 


.08 


<-0.86 


- 


-1.45 + 0.08 


-1.00 + 0.10 


- 


- 


- 


8 


Q1055 + 4611 


3 


,317 


20. 


.34 + 0. 


.10 


- 


- 


- 


>-l,61 


<-1.19 


>-1.97 


<-2.18 


4 


HE1104-1805 


1 


.662 


20. 


.85 + 0. 


.01 


-1.04 + 0.01 


-1.47 + 0.01 


-1.58 + 0.02 


-1.03 + 0.02 


- 


-0.79 + 0.20 


-1.74 + 0.16 


15 


BRIllOS-07 


3 


.608 


20. 


.50 + 0. 


.10 


- 


- 


-2.12 + 0.10 


-1.80 + 0.10 


- 


>-2.37 


- 


8 


BR1117-1329 


3 


.350 


20. 


.84 + 0. 


.13 


-1.18 + 0.13 


-1.36 + 0.13 


-1.51 + 0.13 


-1.26 + 0.13 


- 


>-1.25 


<-2.24 


13 


BRI1202-0725 


4 


.383 


20. 


.60 + 0. 


.13 


- 


- 


-2.22 + 0.16 


-1.77 + 0.12 


- 


>-1.84 


<-2.23 


2,4 


Q1210 + 17 


1 


.892 


20. 


.60 + 0. 


.10 


-0.90 + 0.10 


-1.00 + 0.10 


-1.15 + 0.12 


-0.87 + 0.10 


- 


- 


- 


8 


GC1215 + 3322 


1 


.999 


20. 


.95 + 0. 


.07 


.-1.29 + 0.08 


-1.52 + 0.10 


-1.70 + 0.08 


-1.48 + 0.07 


- 


>-2.56 


- 


8 


Q1223 + 17 


2 


.466 


21. 


.50 + 0. 


.10 


-1.62 + 0.10 


-1.68 + 0.10 


-1.84 + 0.10 


-1.59 + 0.10 


_ 


>-2.76 


-2.60+0.18 


8,10 


Q1232 + 0815 


2 


.338 


20. 


.80 + 0. 


.10 


_ 


_ 


-1.59 + 0.13 


-1.18 + 0.13 


-1.17+0.14 


_ 


-2.10+0.13 


19 


MC1331 + 170 


1 


.776 


21. 


.18 + 0. 


.04 


-1.30 + 0.05 


-1,97 + 0,04 


-2.06 + 0.05 


-1.45 + 0.10 


- 


- 


-1.82+0.10 


8,10 


BRI1346-0322 


3 


.736 


20. 


.72 + 0. 


.10 


- 


- 


<-2.09 


-2.33 + 0.10 


- 


>-2.44 


- 


8 


Q 1351 + 318 


1 


.149 


20. 


.23 + 0. 


.10 


-0.38 + 0.17 


-0.94 + 0.17 


-0,99 + 0.12 


-0.56 + 0.17 


- 


- 


- 


16 


Q1354 + 258 


1 


.420 


21. 


.54 + 0. 


.06 


-1.63±0.16 


-1.82±0.09 


-2.02 + 0.08 


-1.74 + 0.15 


- 


— 


- 


16 


Q1409 + 095 


2 


.456 


20. 


.54 + 0. 


.04 


- 


— 


-2.30 + 0.02 


-2.01+0.02 


- 


-2.13+0.06 


<-3.28 


12 


Q 1425 + 6039 


2 


.827 


20. 


.30 + 0. 


.04 


- 


- 


-1.33 + 0.04 


>-1.03 


- 


— 


-1.53+0.04 


10 


GB1759+7539 


2 


.625 


20. 


.80 + 0. 


.10 


>-1.78 


-1.26 + 0.10 


-1.21 + 0.10 


-0.82 + 0.10 


-0.76+0.10 


>-1.28 


-1.53+0.02 


8,10 


0^1946 + 7658 


2 


.843 


20. 


.27 + 0. 


.06 


- 


- 


-2.53 + 0.06 


-2.23 + 0.06 


<-1.98 


-2.19+0.06 


-3.61+0.07 


8,10 


C^2206-19 


1 


.920 


20. 


.65 + 0. 


.07 


-0.41 + 0.07 


-0.68 + 0.07 


-0.86 + 0.07 


-0.42 + 0.07 


- 






8 


Q2206-199N 


2 


.076 


20. 


.43 + 0. 


.06 


<-1.90 


<-2.19 


-2.61+0.06 


-2.31+0.07 




-1.94+0.09 


<-3.49 


8,12 


C^2212-1626 


3 


.662 


20. 


.20 + 0. 


.08 






<-1.77 


-1.91+0.08 




>-2.18 


<-2.55 


2,4 


LBQS2230 + 0322 


1 


.864 


20. 


.85 + 0. 


.08 


-0.72 + 0.10 


-1.12 + 0.10 


-1.17 + 0.09 


-0.75 + 0.10 








8 


Q2231-0015 


2 


.066 


20. 


.56 + 0. 


.10 


-0.88 + 0.10 


-1.06 + 0.10 


-1.40 + 0.10 


-0.87 + 0.10 








8 


Q 223 7- 0608 


4 


,080 


20. 


.52 + 0. 


.11 






-2.17 + 0.16 


-1.81+0.11 






<-2.17 


2,4 


HE2243-6031 


2 


.330 


20. 


.67 + 0. 


.02 


-1.12 + 0.05 


-1.12 + 0.09 


-1.25 + 0.02 


-0.87 + 0.04 


-0.85 + 0.03 


-0.68 + 0.21 


-1.72 + 0.02 


11 


B2314-409 


1 


.874 


20. 


.10 + 0. 


.20 


<-1.21 


<-1.62 


-1.87 + 0.24 


-1.87 + 0.22 


<-1.66 


>-2.05 + 0.22 




9 


B2314-409 


1 


.857 


20. 


.90 + 0. 


.10 


-1.04 + 0.14 


-1.21±0.13 


-1.32 + 0.14 


-1.05 + 0.14 


-1.00 + 0.18 






9 


Q2343 + 1230 


2 


.431 


20. 


.35 + 0. 


.05 


-0.57+0.06 




-1.19 + 0.06 


-0.65 + 0.06 


-0.83 + 0.08 


>-1.73 


-1.60 + 0.10 


4,17 


Q 2344 + 124 


2 


.538 


20. 


.36 + 0. 


.10 






-1,83 + 0.10 


-1.74 + 0.10 


<-1.36 


>-2.07 


-2.51+0.10 


8,10 


Q2348-01 


2 


.615 


21. 


.30 + 0. 


.10 


<-2.10 


-2.30 + 0.12 


-2,23 + 0.13 


-1.97 + 0.12 








8 


Q 2348- 01 


2 


.426 


20. 


.50 + 0. 


.10 




<-1.46 


-1,39 + 0.10 


-0.69 + 0.10 








8 


Q2348-14 


2 


.279 


20. 


.56 + 0. 


.07 






-2,24 + 0.10 


-1,92 + 0.10 


-2.03+0.10 




<-3.26 


8,10 


Q2359-02 


2 


.095 


20. 


.70 + 0. 


.10 


-0.77+0.10 


-1.55 + 0.10 


-1.65 + 0.10 


-0.78 + 0.10 








8 


Q2359-02 


2 


.154 


20. 


.30 + 0. 


.10 


<-1.07 


-1.37+0.10 


-1.88 + 0.10 


-1.58 + 0.10 
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17) Dessauges-Zavadsky et al. 2002; 18) Petitjean et al. 2002; 19) Centurion et al. 2003. 



are flatter than the oxygen ones in spiral disks as being 
due to differences between the production of O in SNc II, 
which grows with the progenitor mass from 10 Mq to 50 
Mq, and that of S, which also occurs in high mass stars, 
but in a more restricted mass range though. This interpre- 
tation is consistent with chemical evolution models (Timmes 
et al. 1998) using the Woosley & Weaver (1995) yields for 
SNell, that exhibit [S/Fe] strongly decreasing with [Fe/H] 
for [Fc/H] <-2. A third explanation (Izotov & Thuan 1999; 
Izotov, Schracrcr & Charbonnel 2001) for the variations of 
S/Si is that ionization corrections for these elements are im- 
portant, since S and Si are singly ionized and could also oc- 
cur in HII regions (Viegas & Gruenwald 1991); in contrast, 
O is mostly neutral in H I gas. All these caveats should be 
kept in mind when using S as a proxy of a elements. Note, 
however, that the [0/Si] ratio in the sample has suprasolax 
values, [O/Si] median = 0.18, which is well explained by the 



dust depletion scenario, since O is essentially undepleted in 
the ISM. Unfortunately, in most cases the abundance of O is 
very uncertain because it is estimated from saturated lines; 
in this case, we use S as a proxy of O. 



2.2 Dust Correction 

The Si and Fe abundances are particularly important in the 

study of the [a/Fc] pattern, since Si and Fc arc the ele- 
ments tracing the a and iron peak group elements, respec- 
tively, with the largest number of abundance measurements 
for DLAs. However, since Fe and Si are refractory elements, 
they are expected to be depleted into grains. Ideally, one 
would like to gauge the level of dust depletion, and to ap- 
ply dust corrections to the observed abundances in order to 
recover the intrinsic ISM (dust and gas phases) abundances. 
In DLA studies, the classical indicator of dust deple- 
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[Zn/H] 



Figure 1. [Fe/Zn] and [Cr/Zn] vs. [Zn/H] observed in DLAs. The 
dotted line represents the solar value. 




[Sl/H] 



Figure 2. [S/Si] vs. [Si/H] (filled circles) and [O/Si] vs. [Si/H] 
(crosses) observed in DLAs. The dotted line represents the solar 
value. The heavy solid lines represent the robust statistical fits 
to [S/Si] vs. [Si/H]: MM (middle line), USMM (upper line) and 
LSMM (lower Une). 



tion has been the Zn/Fe (or Zn/Cr) ratio. Pettini at al. 
(1994) demonstrated that the majority of DLAs exhibit an 
enhancement of Zn relative to Cr. This overabundance have 
interpreted by them as an indication of dust depletion be- 
cause Zn tends to track Fe in stars of all metallicity (Sne- 
den, Gratton, & Crocker 1991), whereas Zn/Fe is enhanced 
in dust depleted regions of the ISM because Zn is mildly 



refractory while Fe is heavily depleted (Savage & Sembach 
1996, hereafter SS96) 

However, recent stellar analyses have re-examined the 
Zn/Fe ratio in metal-poor and thick disk stars and have 
found enhanced Zn/Fc ratios (Johnson 1999, Prochaska ct 
al. 2000, Umeda & Nomoto 2002). These surveys find an 
enhancement of ~ +0.1 in the Galactic thick disk stars 
([Fe/H] « —0.5) (Prochaska et al. 2000) and a mean en- 
hancement of « 4-0.2 to -f-0.3 in metal-poor stars ([Fc/H] < 
—2) (Johnson 1999). These results indicate the possibility of 
a nucleosynthetic contribution to the Zn/Fe enhancement of 
a given DLA. On the other hand, in the Galactic ISM, the 
suprasolar [Zn/Fc] and [Si/Fc] ratios provide evidence for 
dust depletion of Fe and Si , but the dust depletion pat- 
tern depends on the type of environment through which the 
sight line passes (SS96). SS96 consider four different dust 
depletion patterns corresponding to four types of Galactic 
ISM environment: (1) cool clouds in the Galactic disk (CD), 
(2) warm clouds in the disk (WD), (3) disk plus halo clouds 
(WHD), and (4) warm halo clouds (WH). 

Applying dust corrections is a very complex issue, and, 
therefore, we considered a variety of dust correction mod- 
els, based on the work of Vladilo (2002a) . Following Vladilo 
(2002b), we apply the dust correction models of Vladilo 
(2002a) to our sample of DLAs with [Zn/Fe], [Fe/H] and 
[Si/H] measurements, in order to obtain the dust-corrected 
[Fe/H] and [Si/Fe] abundances. Our prescriptions, however 
differ from those of Vladilo (2002b) in that: i) we allow for 
a larger scatter in the nucleosynthetic contributions to the 
Zn/Fe enhancement by assuming a range for the intrinsic 
[Zn/Fe] ratio of [Zn/Fe]i = 0.0-0.2, where Vladilo (2002b) 
considers [Zn/Fe]i = 0.0 — 0.1; ii) we take as reference me- 
dia the WD and the WH, whereas Vladilo (2002b) consid- 
ers the WD and the CD. In fact, it is not appropriate to 
consider the cool clouds as a reference for the DLA environ- 
ment, because the cool clouds exhibit levels of dust depletion 
([Fe/Zn], [Cr/Zn], [Si/Zn]) much higher then those observed 
in DLAs, and the molecular hydrogen fractions are typically 
low in DLAs, in contreist to the large amount of molecules 
in the cool clouds. In addition, Savaglio, Panagia & Stiavelli 
(2000) have found, for 37 DLAs, that the majority (~ 60%) 
show a WH-likc depletion pattern, then ~ 30% and ~ 10% 
have a WHD-like and WD-like depletion pattern, respec- 
tively. Finally, Hou et al. (2001) have found evidence for a 
dust depletion, for DLAs with large metal column density 
([Zn/H]+logN(H I)> 19.8), at a level (~ 0.4 dex) compara- 
ble to that observed in the WH clouds. 

We summarize below the dust correction method of 
Vladilo (2002a). The dust correction method aims at ob- 
taining the depletion of a given element X in the DLA: 



6x = [X/HUs - [X/H]i 



(1) 



where [X/H]obs is the observed (gas) abundance and [X/H]i 
the intrinsic ISM abundance of the DLA. The depletion 5{X) 
is related to the fraction in dust of the element X, fx = 
Nx,d/Nx,i, through: 



fx =1-10" 



(2) 



where Nx,d is the column density of atoms X in the dust 
and Nx,i the column density in the ISM (gas plus dust) of 
the DLA. 

The fraction fx is the relevant physical quantity, which 



© 2003 RAS, MNRAS 000, 000-000 



Nature of DLAs and chemical evolution 7 



should scale with the amount of dust and with the abun- 
dance of X in dust. To derive this relation we express all 
abundances by number and relative to the abundance of a 
reference element Y (here we adopt iron as Y). Therefore, 
we define the relative abundance of X in the ISM as: 

ax = Nx,i/NY,i , (3) 

which is related to [X/Y]i through ax = lO'^/'^'-. Wo also 
define the reference dust fraction (or dust-to-metals ratio as 
in Vladilo 2002a) in the ISM as fy- 

We further assume that fx depends on ax and fy as: 



(4) 



in which {l + r)x) and {ex — 1) are the power law indices of 
the power laws in fy and ax , respectively. Eq. (4) allows us 
to obtain a general expression for fx : 



fx 



fy 



fy 



ax 

ax,m 



(5) 



If fx,m is known for some reference medium m, and, if rjx 
and ex are known or assumed, we can obtain fx- Here we 
use as reference environment m the warm disk (WD) of the 
Galaxy. Our values of fx,m correspond to 5x,m ~ —0.43, 
-1.095, and -1.215, for Si, Cr and Fe, respectively (Sombach 
& Savage 1996), and 6x,m = -0.19 for Zn (Roth & Blades 
1995). 

In Vladilo (1998), fx was assumed to scale linearly with 
fy. However, this is correct when the dust composition does 
not vary with fy , a situation described by r/x = 0. However, 
studies of the Galactic ISM indicate that the dust compo- 
sition changes with fy. In particular, for elements that are 
less depleted than the iron (e.g. silicon) the fraction in dust, 
in general, decreases faster than that of the iron for less 
depleted media, implying r/x > 1. 

On the other hand, ex = implies that the dust com- 
position is insensitive to variations of metallicity. However, 
one expects that the dust fraction would increase monoton- 
ically with the metallicity. In a simple physical picture, the 
two-body process of adhesion of an element onto dust grains 
should be proportional to the product of the densities of that 
element and the dust grains. If there are fewer atoms avail- 
able, then the adhesion rate will bo correspondingly lower, 
ex = 1 follows from this picture, implying that the dust 
fraction scales linearly with the metallicity of the medium. 
Values of ex > 1 could occur if three-body processes are 
important in the chemistry of dust production, but this is 
never the case in the low densities of the ISM. In addition, 
ex > 1 not supported by studies of the ISM of the SMC 
(Vladilo 2002a). Therefore, the physically meaningful range 
for ex is < ex < 1. 

The value of rjx can be obtained from the compari- 
son of knowledge of the dust depletion pattern in the local 
Galactic ISM, considering two different environments m and 
n. We assume that the abundances are constant in the so- 
lar neighbourhood. This implies that ax,m = ax,n and, as 
a consequence, the scaling law, eq. 5, can be written in a 
simplified form: 



fy 



fY,n 



{1+'7X) 



(6) 



media used in the determination of rjx. For the warm halo, 
Sembach & Savage (1996) derive 5x,wh = —0.26, -0.51, 
and -0.64, for Si, Cr and Fe, respectively. In their study of 
Zn depletion in the Galactic ISM, Roth & Blades (1995) 
included one sightlinc of Sembach & Savage in the WH cat- 
egory, that toward HD167756, and, in addition, the WDH 
class sightline toward HD149881. They find for the sight- 
lines toward HD167756 and HD149881, 8zn,WH = -0.05 
and -0.11, respectively. Accordingly, in this work we assume 
5zn,WH ~ —0.10. However, assuming 5zn,WH = —0.05 has 
little effect on our results (see Table 2). The corresponding 
values of jjx are 0.685, 0.450 and 1.760, for Si, Cr and Zn, 
respectively (=4.922 for Zn if 5zn,WH = -0.05). 

Once wc have determined or estimated ex and rjx, 
we can proceed in our attempt of dust correction of the 
abundances, by considering the relation between the ob- 
served [X/H]o6s and the intrinsic abundance [X/H]i (or 
ax = iqI'^/^I' in the equations below), which can be de- 
rived from Eq.s (1) and (5) above: 



x 




x- 


.1. 


ohs 


H. 



-|-log{l-p 



/X.r.} , (7) 



where p — fy / fy,m is the normalized (to fy,m of the refer- 
ence medium m) dust-to-metals ratio. 

In order to recover the intrinsic abundances we first 
need to determine p. With this task, we apply Eq. (7) to the 
generic element X and to the reference element Y to obtain 
the relation: 



fx,m ax 
fy.m lol^l 







(8) 



If we have a measurement of [X/Y]oi,s and a guess on 
[X/Y]i, Eq. (8) can be solved to find the root p. 

Once p is known, wc can determine the intrinsic relative 
abundance [X/Y]i of a element X different from that used 
in Eq. (8) from the expression 



ax 



■ p^'^"^' fx,^ ax + (p/y.™ - 1) loI^Li.. = 



(9) 



Here, the warm disk and the warm halo are the two 



which is just a rearrangement of Eq. (8), with ax as the 
unknown variable instead of p. 

In applying the dust correction procedure of elemental 
abundances described above, we follow Vladilo (2002a), and 
consider eight dust models, based on our guesses on r/x and 
[Zn/Fe]i. We consider [Zn/Fe]i as an input parameter to 
which we assign values 0.0 or -1-0.2 dex, in order to take into 
account a possible nucleosynthetic Zn/Fe enhancement. We 
adopt Tjx = or 77x = 1 as two limiting possibilities for r)x . 
As rjx is needed both for zinc in Eq. (8) in order to obtain 
p, and for the generic element X in Eq. (9) , we consider that 
possibility that rjzn and rjx arc different. We have, therefore, 
eight possibilities: SOO, SIO, SOI, Sll, EOO, ElO, EOl, and 
Ell, where the letter S or E denotes solar [Zn/Fe]i = 0.0 
(S) or Zn enhanced [Zn/Fe]i = 0.2 (E). The number code 
indicates the values of the pair {rjzn,rjx) = (0.0,0.0) (00), 
(1.0,0.0) (10), (0.0,1.0) (01), and (1.0,1.0) (11). We should 
note that SOO and SIO {= SO) and SOI and Sll {= SI) are 
indistinguishable. In the case of iron abundances, SO and 
SI arc indistinguishable, and also EOO and ElO, and EOl 
and Ell. For each dust correction model, the correction is 
applied to Si, Cr, Fe and Zn in the systems for which the 
observed [Zn/Fe] is available. 
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Figure 3. Upper panel: Dust-corrected [Si/Fe] vs. [Fe/H] for 
DLAs (filled circles). Also show [S/Fe] vs. [Fe/H] data (open cir- 
cles). Lower Panel: [N/a] vs. [a/H] for DLAs and dwarf galaxies. 
The dwarf galaxy data arc [N/O] vs. [O/H] from van Zee et al. 
(1997a, b) (crosses) and Izotov & Thuan (1999) (plus signs). For 
DLAs, wc take as a proxy for a elements O or S, or, when not 
available, dust-corrected Si. The dotted horizontal lines represent 
the solar value, the dotted vertical line in the bottom panel the 
value [a/H]=-1.3 (see text), and the heavy solid lines are the sta- 
tistical fits to [Si/Fe] vs. [Fe/H] and [N/a] vs. [a/B] of DLAs. The 
bottom panel also shows the lower bound envelope to the [N/a] 
vs. [q/H] distribution of DLAs (dashed line). 

Table 2 shows, for each system, the observed abundance 
ratios [Zn/Fe] and [Si/Fe] and the observed abundance of 
the iron [Fe/H], the corrected abundances [Fe/H] cor, the 
corrected abundance? ratios [Si/Fc]cor, for each dust cor- 
rection model, and the averaged corrected [Fe/H] cor and 
[Si/Fe] cor-, ([Fe/H]) and ([Si/Fe]) (the average has been done 
over the eight possibilities SOO, SOI, SIO, Sll, EOO, EOl, 
ElO and Ell). The error bars of ([Fe/H]) and {[Si/Fe]) en- 
compass the range of average values predicted by each dust 
model and their respective la errors. All models assume 
[Zn/H]vKif = —0.1 for the Galactic Warm Halo medium, 
but adopting [Zn/H]vF_r/ — —0.05 changes very little our re- 
sults even for the system with the largest [Zn/Fe] obs, the 
Zabs = 2.095 DLA toward Q2359-02 (see Table 2). 



2.2.1 The [Si/Fe] ratio 

The trend of the dust-corrected [Si/Fe] ratio with mctallic- 
ity {[Fe/ H]cor) is analysed with robust statistical methods 
(curves MM, USMM and LSMM in Figure 3a). In our sam- 
ple, before any dust correction, the [Si/Fe] ratio is typically 
suprasolar, with [Si/Fe]median = 0.32. In the corrected sam- 
ple, the [Si/Fe] ratio is lowered to [Si/Fe]medmn = 0.16, but 
remains suprasolar. Notice that the original dust correction 
procedure of Vladilo (1998) is equivalent to the case SO, 
which in general gives a larger downward correction for the 



[Si/H] ratio. After we perform a careful dust correction, we 
conclude that DLAs do exhibit [Si/Fe] enhancement, and 
there is a trend of slight increase of [Si/H] with decreas- 
ing metallicity. resembling that found for metal-poor stars 
of the Galaxy. This suggests that a number of DLAs may 
have had a chemical evolution similar to our Galajcy, sup- 
porting the scenario where these systems axe the progen- 
itors of disks. However, there are DLAs with nearly solar 
[Si/Fe] ratio even at low metallicities, for which the disk 
picture does not apply. This case implies a slow star forma- 
tion rate, which allows the SNe la to enrich the ISM with 
iron peak group elements. This behaviour is similar to that 
predicted by chemical evolution models for dwarf galaxies 
with bursts of star formation separated by quiescent peri- 
ods (Matteucci et al. 1997), which allows that class of DLAs 
to be explained as (proto)dwajrf galaxies. In summary, the 
[Si/Fe] vs. [Fe/H] distribution of DLAs does not exclude 
disks as DLA sites, but, if disks are present among DLA 
hosts, implies that they represent one component within a 
multi-population scenario. 

2.2.2 The [N/a] ratio 

The [N/a] ratio also allows one to constrain star formation 
times scales because of the delay of the injection of N with 
respect to a elements into the ISM. A nearly solar ratio 
([N/a] = 0) would indicate a slow star formation, during 
which there is enough time for production of N. Figure 3b 
show the trend [N/a] vs. [a/H] for DLAs and dwarf galaxies. 
The spread in the values of [N/a] observed in DLAs is larger 
than that observed in our Galaxy for stars in the same range 
of metallicity and in dwarf galaxies, and can be attributed to 
different formation epochs or chemical evolution histories for 
these systems. The high metallicity ([a/H] > —1.3) DLAs 
exhibit values similar to those of dwarf galaxies giving sup- 
port to the suggestion that DLAs could be the progenitors 
of dwarf galaxies. At lower metallicities, however, the [N/a] 
decreases, and there are DLAs with [N/a] values below that 
of any dwarf galaxy, which could be very young objects, for 
which there was no time for injection of N into the ISM. 

Recently, Prochaska et al. (2002b) have claimed that 
the observed [N/a] ratios in DLAs present a bimodal distri- 
bution in the [N/a] vs. [a/H] plot, with a small sub-sample 
exhibiting significantly lower N/a ([N/a]~ —1.5) than the 
majority of systems. The authours consider that a way of 
naturally reproducing a bimodal distribution of [N/a] val- 
ues is a top heavy IMF, thus reducing the nitrogen con- 
tribution of intermediate mass stars (IMS). The bimodal- 
ity would severely challenge other explanations: the low - 
DLAs [N/a] are observed during a transient period occur- 
ring within 250 Myr of a star burst and before the IMS have 
released their nitrogen into the ISM; N yields of IMS sig- 
nificantly decline with decreasing metallicity; the low [N/a] 
reflects the nucleosynthesis of Population HI stars. Notice 
that a test of the top heavy IMF explanation would be that 
it predicts a high [a/Fe] ratio. 

However, in our sample, the bimodality in the [N/a] vs. 
[a/H] plot is less apparent. Notice that we consider O, S, and 
Sicor as the a-elements, while Prochaska et al. (2002b) use Si 
as a proxy of a-elements and do not consider dust-correction. 
In our sample, there are three systems with determinations 
of N and a abundances, which have low (~ —1.5) [N/a]. 
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Table 2. Corrected abundance of Fe and Si. 
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" The dust models assume [Zn/H]wif = —0.05 (see text). 



Only one system, the DLA at z =2.843 toward Q 1946+7658 
which happens to be that with the lowest metallicity in Fig- 
ure 3, is isolated from the rest of the systems. The other two 
low [N/a] DLAs seem to be contiguous to the DLAs with 
higher [N/q]. 

Independently from a segregation of a low [N/a] sub- 
sample, to be confirmed by further observations, there seems 
to be a lower bound envelope in the [N/a] vs. [a/H] plot, 
extending in an arc from ([a/H], [N/a]) ~ (—1.5,-1.5) to 

~ (-0.7, -0.9). The two low [N/a] DLAs with [a/H] 1.4 

are located along this envelope, forming a sequence with 



higher [N/a] DLAs. This feature of the [N/a] vs. [a/H] dis- 
tribution is interesting from the point of view of the models 
discussed in this work, since it corresponds to a ma^ximum of 
the specific star formation rate v (i.e. the inverse of the star 
formation time scale), which is a central ingredient of the 
chemical evolution process. Also the upper limits on [N/a] 
are consistent both with non-bimodality and an upper limit 
for v. 
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log N(HI) + [Si/H] or log N(HI) + [Zn/H] 

Figure 4. The trend of [Zn/Fe] with the metal column density 
F(X) for X= Zn, Si, and Sicor- The dotted line represents the 
solar value. The thick solid lines are the statistical fits to the 
DLA data. 
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Figure 5. The trend of [Zn/Fc] with [Si/H] for the DLA sample 
with Zn, Fe and Si measurements. The dotted line represents the 
solar value. The thick solid lines are the statistical fits to the DLA 
data. 



2.3 Estimating the dust depletion for DLAs 
without [Zn/Fe] ratios 

As any nucleosynthetic contribution to the [Zn/Fe] ratio is 
small, not reasonably exceeding 0.3 dex, the [Zn/Fe] en- 
hancement is a reliable indicator of dust depletion. This is 
confirmed in Figure 4 from the good correlation between 
the observed [Zn/Fc] ratio and the metal column density 
F(Zn)=log N(HI)-|-[Zn/H], which is used in the definition of 
"dust filter" by Pramtzos & Boissier (2000) , based on Boisse 
et al. (1998). The trend of increasing [Zn/Fe] with F(Zn) 



seems to fiat beyond F(Zn) somewhat below F(Zn)=20.0. 
One could wonder if obscuration effects due the higher col- 
umn density of metals would not prevent the highest de- 
pleted media from been seen. It is interesting to note that 
this apparent break in the trend occurs essentially at the 
same value of F(Zn)=19.8, beyond which Hou et al. (2001) 
find evidence for dust depletion in DLAs. 

If we consider the equivalent definition of F(Zn) for the 
silicon, F(Si)=log N(HI)-|-[Si/H], the correlation is weaker. 
The reason for that, is that the weakly depleted zinc allows a 
much more reliable determination of the total metal column 
density than the Si, which is mildly depleted. However, if 
wc use the dust corrected Si for the individul DLAs, we 
obtain a much more clear trend of [Zn/Fc] with F(Sicor). The 
correlation [Zn/Fe]-F(Sicor) allow us to be more confident 
about the dust correction method used in this paper. In 
addition, the trend of increasing [Zn/Fe] with F(Sicor) flats 
for F(Sicor)<^ 20, which is the same behaviour exhibited by 
F(Zn). However, in this case, the break seems to occur at 
a value of F(Sicor) that is slightly higher that value of the 
break of F(Zn). This would be expected if there is a modest 
enhancement of the [Si/Fe] ratio. 

The considerations above allow us to try to extend the 
sample of DLAs with estimated corrected abundances. In 
fact, only w 40% of the DLAs in our database have been se- 
lected to be applied the dust correction described in the last 
subsection, and important pieces of information could be 
missed. Wc would like to estimate intrinsic (dust-corrected) 
[Fe/Si] ratios for systems without Zn measurements. In the 
following, we will try to use silicon as an indicator of dust 
depletion, or, equivalently, of the Zn/Fe ratio. Although sil- 
icon is a refractory element, its depletion is mild in lightly 
depleted regions of the ISM, and we have found that the cor- 
rections to be applied to Si were never large. The effects of 
dust depletion may hide significant trends with metallicity, 
and these effects can be reduced if we use a little refractory 
metallicity indicator. In fact, dust depletion has spoiled the 
correlation between [Zn/Fo] and F(Si), as wc have already 
seen from the less clear trend of [Zn/Fe] with F(Si), when 
we use the observed (and therefore affected by dust) [Si/H] 
instead of the dust-corrected [Si/H]. However, as wc can see 
from Figure 5, there is a clear correlation between [Zn/Fe] 
and [Si/H], oven using the observed [Si/H]. 

As the [Si/H] metallicity has been revealed to be a good 
predictor of the [Zn/Fc] ratio, we will extend our sample of 
DLA with estimated dust depletion to include those DLAs 
which have a value of F(Zn) too small for the Zn lines to 
be detected. Hereafter, this sample will be called the "no 
Zn sample". It would allow us to explore a domain with 
smaller metallicitics on average than the typical [Fe/H]^ 
—2.0 region defined by the original sample with Zn and Fe 
to which we applied the dust correction procedure. 

The [Zn/Fe] ratio, however, at least up to a metallic- 
ity [Si/H]= —0.8, varies between zero and values close to 
0.8. Below [Si/H]=-1.0, the dispersion of [Zn/Fe] becomes 
smaller and the largest values are bounded by the USMM of 
the statistical trends, with the exception of the Zabs = 1.776 
DLA toward MC1331-I-170. This DLA, however, may be pe- 
culiar, since, in the sample used by Prochaska (2002) to 
investigate abundance variations within DLAs, it presents 
the second largest deviation from uniformity. In order to 
obtain the no Zn sample, we assume, for the systems with 
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Figure 6. Upper panel: The [Si/Fe] vs. [Fe/H] distribution for 
the no Zn sample. Lower panel: Comparison of the Si/Fe] vs. 
[Fe/H] distributions of the no Zn sample (represented by poly- 
gon enclosing this sample) with the sample with Zn, Fe and Si 
measurements. 

measurements of Fe and Si, but not of Zn, and satisfying 
[Si/H]< —1.0, that the [Zn/Fe] ratio varies between the 
USMM and the LSMM of the statistical trends shown in 
Figure 5, and then apply the dust correction method de- 
scribed in Section 2.2. We also use some useful upper limits 
whenever they define a narrower range of [Zn/Fe] than that 
delimited by the USMM and the LSMM. It is the case of the 
DLAs toward PKS0528-2505 at z =2.141 ([Zn/Fe]< 0.17), 
toward Q0836-I-11 at z =2.465 ([Zn/Fe]< 0.27), and toward 
Q2348-01 at z =2.615 ([Zn/Fe]< 0.13). The resulting esti- 
mates of the dust corrected abundances are shown in Fig- 
ure 6. In Figure 6, the polygon encloses the no Zn sam- 
ple. In drawing the polygon, we have excluded the DLA at 
z =1.874 toward B2314-409, which has an unusually low 
observed [Si/Fe] ratio (Elhson & Lopez 2001, 2002). Due 
to the uncertainties in the dust correction method and the 
two-step nature of the process used in building this sample, 
assigning dust-corrected abundances has no sense for indi- 
vidual DLAs but only for the sample as a whole. The no Zn 
sample seems to be an extension toward low metallicities of 
the sample with measurements of Zn, Fe and Si. It is not 
possible to know if the trend of increasing [Si/Fe] with de- 
creasing [Fe/H] continues at [Fe/H]< 2 or if [Si/Fe] reaches 
a plateau. For this sample, [Si/Fej^edian = 0.204, an en- 
hancement only ~ 0.05 dex higher than that of the sample 
with Zn, Fe and Si measurements. 



3 THE CHEMICAL EVOLUTION MODELS 

In order to study the evolution of the metal content in DLAs, 
we use two classes of models: one-zone chemical evolution 
models and chemodynamical models. The one-zone chem- 



ical evolution models are used to test the disk or dwarf 
galaxy scenario for DLAs, whereas the chemodynamical 
model checks for the dwarf galaxy scenario. 

For each model class, a star formation law is specified. 
The stars formed follow a Salpeter IMF between 0.1 and 100 
Mq. Chemical evolution occurs as the stars formed out of 
the ISM evolve and eject gas back into the ISM. The stars 
are assumed to die either as supernovae (SNe) or as plan- 
etary nebulae, when instantaneous ejection of mass, metals 
and energy occurs. The evolution of the abundances of He, 
C, N, O, Mg, Si, S and Fe is calculated by solving the basic 
equations of chemical evolution (see Matteucci & Tornambe 
1987). We do not assume instantaneous recycling approx- 
imation for the chemical enrichment, but instead we take 
into account the delays for gas restoring from the stars due 
to the main-sequence lifetimes. Instantaneous mixing with 
the ISM is assumed for the stellar ejecta.. The models start 
with an entirely gaseous protogalaxy or infalling gas with 
primordial chemical abundances {Y = 0.24, Z = 0). 

In this work, we use metallicity dependent yields for SNe 
II, SNe la, and IMS. For SNe II, we use the metallicity de- 
pendent yields of Woosley & Weaver (1995), which are given 
for stars of mass M = 12, 13, 15, 18, 20, 22, 25, 30, 35 and 
40 M0 and metallicities Z/Zq = 0, 10"*, 10"^ 10"^ and 
1. Yields of SNIa resulting from Chandrasekhar mass white 
dwarfs, are taken from Iwamoto et al. (1999): their models 
W7 (progenitor star of initial metallicity Z=Zq) and WTO 
(initial metallicity Z=0). The yields for IMS (0.8 - 8 Mq), 
with initial Z=0.001, 0.004, 0.008, 0.02 and 0.4, are from 
van den Hoek & Groenewegen (1997) (their variable ijagb 
case). For more details of the nucleosynthesis prescriptions, 
see Friaga & Terlevich (1998, hereafter FT98). 

Our choice of the mass (or surface mass density) and 
physical size of the model galaxies was based on the scenario 
in which the amount of gas in DLAs constitute a reservoir 
for the formation of present day galaxies, in view of the co- 
incidence of the amount of mass in DLAs and the baryonic 
stellar mass in galaxies today (Storrie-Lombardi & Wolfe 
2000). Therefore, the column density of the models should 
be not far from the baryonic column density of the present 
day galaxies being modelled. At the same time, the H I col- 
umn densities predicted by all the models satisfy both our 
DLA definition, log N(HI)> 20.0, and the high column den- 
sity cut-off log N(HI)^ 22, observed in DLAs, and explained 
by Schaye (2001) as due the conversion of atomic hydrogen 
to molecular hydrogen in clouds with higher total hydrogen 
(HI-I-H2) column density. Therefore, in comparing the obser- 
vations with the predictions of the models, we will consider 
not only the abundances and abundances ratios, but also 
the H I column densities. 



3.1 The one-zone chemical evolution models 

Two types of one-zone chemical evolution models are used 
to explore the two most popular hypothesis for the nature 
of DLAs: disks (Prochaska & Wolfe 1997), and dwarf galax- 
ies (Centurion et al. 2000, Molaro et al. 2001). One- zone 
models with galactic winds are used for dwarf galaxies, and 
multi-zone chemical evolution models with infall are used to 
represent disks. 
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3.1.1 The dwarf galaxy model 

The dwarf galaxy model is the classical model with galactic 
winds (Matteucci 1994), first proposed by Larson (1974). 
In this model a wind is established at a time tw when the 
thermal energy of the gas becomes higher than the potential 
energy, i?th,sjv > -Eb.gas- The wind instantly sweeps out all 
the gas of the galaxy and the star formation ceases. Later, 
processes of stellar mass loss restore the ISM. 

Most of the models represent a galaxy with initial gas 
mass Mg = 10® M©, but there is also a model with 
Mq ~ 10^ M0. The galaxies are inside a dark halo with 
a mass 7.5 times Mq- The star formation rate (SFR) is pro- 
portional to the present mass of gas and characterised by the 
specific star formation rate v (the inverse of the star forma- 
tion time scale). The fact that these models are simple, pa- 
rameterized only in terms of u makes them useful for gauging 
the effect of the star formation rate in the chemical evolu- 
tion of DLAs. We run models for several values of v ranging 
from 0.1 to 19 Gyr~^. = 19 Gyr~^ follows from assuming 
Ma = 10'^ Mq in the relation = 8.6(Mg/10^^ M0)-°-"^ 
used in Matteucci & Tornambe (1987) and later works. The 
values of < 3 Gyr~^ reflect the typically low SFR of DLAs, 
as inferred from observations. 

Interesting limits on the SFR of DLAs are given by 
imaging observations of the Zab.i = 3.386 DLA toward 
Q0201-I-1120 (Ellison et al. 2001a). There are only two faints 
sources near the quasar with photometric redshifts consis- 
tent with the redshift of the DLA, the LBGs G2 and oM6 
separated from the QSO by 4 and 2.9 arcsec. Spectroscopy 
of oM6 (7?.=23.97 mag) revealed a prominent Lyman a emis- 
sion line at Zem ~ 3.645. Therefore, oM6 is not the DLA. 
Rather, it is associated to the QSO (the difference in redshift 
between the oM6 and Q020H-1120 is only -1-465 km s~^). 
G2, the remaining candidate for DLA, is so faint (7?-=25.3) 
that there is no spectroscopy available. However, assuming 
that G2 is at the redshift of the DLA, we can estimate its 
SFR from its TZ magnitude, because, at a: ~ 3, the TZ fil- 
ter samples the rest-frame far-UV spectrum produced by 
recently formed O and B stars. Assuming that G2 has a 
typical DLA metallicity of [Fe/H]= —1.5 and that is un- 
dergoing continuous star formation for 10® yr, we obtain, 
following Friaga & Terlevich (1999), a SFR of 6.2 Mq yr~\ 
This corresponds to a model with Mq ~ W'' Mq and u — 6 
GyT-^ (or Mg = 2 X lO'' M© and = 3 Gyr'^). However, if 
G2 is not the DLA, then the absorber must be fainter than 
TZ « 26.0, and the corresponding SFR <, 3 Mq yr~^, or 
z/ < 3 Gyr-i for Mq = 10® Mq. 



3.1.2 The disk model 

The disk model belongs to the class of the models for the 
Galaxy with one infall episode of pristine gas. Following 
Matteucci & Frangois (1989), the SFR depends on the galac- 
tocentric radius r through the total mass surface density 
cr{r,t) and the gas fraction G(r,t) as: 



*(r,t) = u 
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The normalisation D of the SFR law is set at the solar po- 
sition (f = 8 kpc) and at a galaxy age to = 13 Gyr (cor- 
responding to a galaxy formation epoch zgf = 10). We 



assume i/ — 0.5 Gyr^^ and xsf = 1/2 in the adopted SFR 
law, in which the specific star formation rate depends on a 
power of the gas density v oc p^'^*'. The infall rate into the 
disk is Eissumed to decrease with time as e~*^'^°, with the 
infall time scale td given by: 

TD {r) = 1 + 7/6(r-2), for 2 < r < 8 kpc 

Td (r) = 8 -h 0.5(r-8), for r > 8 kpc 
In this way, td — 8 and 1 Gyr at the solar position and 
at r = 2 kpc, respectively (Chiappini, Matteucci & Graton 
1997), and td ~ 13 Gyr at r = 18 kpc. An exponential sur- 
face mass density profile is adopted for the disk, a oc e~^^'''^ , 
with ra = 2.6 kpc (Boissier & Prantzos 2000), normalized to 
aifJa) = 50 Mg, pc*^ (Chiappini ct al. 1997). We do not 
consider any threshold for star formation, differently from, 
for instance, Chiappini et al. (1997), who adopt the thresh- 
old Eth = 7 Mq pc~^. Adopting a threshold would impact 
on the outer regions of the disk, where it would induce a 
starburst-Iike star formation regime (Calura et al. 2002). 

The disk model, the column density of HI implied by 
the present day total surface mass density is 4.7 x 10^^ cm~'^ 
at the solar radius and 1.0 x 10^° cm~^ at r = 18 kpc, which 
is consistent with the range of N(HI)values in DLAs. During 
the evolution of the disk (not in the very early phases, i.e. 
t < TD{r), when there is still very little gas in the disk), 
the value of the column density will vary by a factor of few 
around the values above. 



3.2 The chemodynamical model 

We investigate the scenario in which the DLAs are dwarf 
galaxies with the chemodynamical model of FT98, in which 
a single massive dark halo hosts baryonic gas that will fall 
toward the centre of dark halo and will subsequently form 
stars. The code combines a multi-zone chemical evolution 
solver and a ID hydrodynamical code. The system, assumed 
to be spherical, is subdivided in several spherical zones and 
the hydrodynamical evolution of its ISM is calculated. The 
equations of chemical evolution for each zone are then solved 
taking into account the gas flow, and the evolution of the 
chemical abundances is obtained. A total of «100 star gen- 
erations are stored during 13 Gyr for chemical evolution 
calculations. We assume that at a given radius r and and 
the time t, the specific star formation rate v{r,t) follows a 
power-law function of gas density (p): i'{r,t) = v{p/ poY^^ , 
where po is the initial average gas density inside the core ra- 
dius of the dark halo (r^), and v is the normalization of the 
star formation law. We include inhibition of star formation 
for expanding gas (V.m > 0) or when the density is too low, 
and, therefore, the cooling is inefficient (i.e., for a cooling 
time tcoo ~ {3/2)kBT / iJ.mHl^{T)p longer than the dynami- 
cal time tdyn = (37r/16 G p)^/^) by multiplying i> as defined 
above by the inhibition factors {1 + tdyn max(0, V.u))~^ and 
(1 + tcoo/tdyn)~^ ■ A characteristic of these models is that 
several episodes of inflow and outflow occur simultaneously 
ar different radii. The chemodynamical model for spheroids 
was used to investigate the relation between young elliptical 
galaxies and QSO activity (FT98), the absence of passively 
evolving elliptical galaxies in deep surveys (Jimenez et al. 
1999), LBGs (Friaca & Terlevich 1999), Blue Core Spheroids 
(Friaga & Terlevich 2001), and the coupled spheroid and 
black hole formation (Archibald et al. 2002). 

We have extended the model family used in previous 
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works (FT98, Friaga & Terlevich 1999, 2001) toward lower 
masses to study dwarf spheroids. The model dwarf galaxy 
has an initial gas mass Mg = 10^ M©, and the normal- 
ization of the SFR was set j/ = 3 or 1 Gyr~^, rather than 
10 Gyr~ as for luminous olhpticals (in order to reproduce 
the suprasolar [Mg/Fe] ratio of giant ellipticals). The scaling 
laws used in FT98, which are valid for giant or intermedi- 
ate ellipticals, were changed to forms more appropriated for 
dwarf galaxies, using as reference the "bright dwarf ellipti- 
cals" defined in Bender, Burstein & Faber (1993): 1) rn = 1 
kpc was kept as in the 10^° Mo model, since dwarf galaxies 
have a larger radius with respect to their masses than bright 
eUipticals; 2) rudail rh = 14 instead of 28, used for himi- 
nous/intcrrncdiatc clfipticals; 3) the Faber- Jackson relation 
was redefined in order to fit the pair NGC205-NGC3605. 
In FT98, the stellar population of the model galaxy fol- 
lows a Faber-Jackson relation, cr« = 200(LB/ii3)^''* km s~^, 
with Lb related to Mg through [M/Ls] = 10, the mass- 
light ratio typical of an L* galaxy, or, equivalcntly, u* = 
150.5(M/10" UVofl^ km s"^ Here, we take the dynam- 
ical masses and velocity dispersions determined by Ben- 
der at al. (1993) for NGC205 and NGC3605 to obtain 
cr. = 188.3(M/10" M0)°-^'^^ km s'^ The initial HI column 
density of the model is N(HI)= 1.65 x 10^^ cm"^ through 
the centre of the galaxy. 



4 INVESTIGATION OF DLAS WITH 
CHEMICAL EVOLUTION MODELS 

4.1 One-zone model for dwarf galaxies 

The predictions for [N/a] of the classic wind model for dwarf 
galaxies are compared to the observed ratios both in DLAs 
and dwarf galaxies in Figure 7. The observed [N/a] in DLAs 
shows a larger dispersion than that found for Galactic stars 
in the same range of metallicity and for dwarf galaxies. The 
different values of [N/a] for a given [a/H] can be explained 
by the different time scales for production of N and a ele- 
ments, as one cam see in the sequence of models with star 
formation ranging from almost quiescent {y = 0.1 Gyr~^) to 
violent (;/ = 19 Gyr^^). Most dwarf galaxies seem to require 
a specific star formation v between 0.3 and 3 Gyr~^, while 
the typical v for DLAs seems to be somewhat lower, with a 
significant number of DLAs with v « 0.1 Gyr~^. For both 
populations, v ftZ Gyr"'^ seems to be an upper limit for the 
efficiency of the SFR. 

In our models, when the galactic wind occurs, the metal- 
licity of the system is too high for a DLA. The time of the 
onset of the is galactic wind is = 0.16, 1.47 and 5.66 Gyr, 
for 1/ = 19, 3 and 1 Gyr^^. The models with v = 0.3 and 0.1 
Gyr~^ do not develop galactic wind until t = 13 Gyr. The 
effect of the galactic wind is illustrated by the model with 
Mg = lO'^ Mq. When the galactic wind occurs at 0.3 Gyr, 
the star formation stops, and, consequently, the O release 
by SNll; however, the IMS keep ejecting N to the ISM. This 
is the reason for the upturn of the Mq ~ 10^ M© model 
evolutionary track at [a/H]w 0.1, away from the Mg = 10^ 
Mq model path. 

The predictions for [Si/Fc] of dwarf galaxy models with 
different values of v are compared to the observed values in 
Figure 8. Most systems have [Si/Fe] ratios below the predic- 
tions of any model. The classic wind model fails in explaining 
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Figure 7. [N/a] vs. [ck/H] observed in DLAs and dwarf galaxies 
compared with one-zone model predictions for dwarf galaxies. The 
solid lines represent the grid of models with v = 0.1, 0.3, 1, 3, 19 
Gyr~^ (left to right), and the dotted line (barely seen) a model 
with V = 3.0 Gyr~^ and M = 10^ M©. Several evolutionary 
times arc indicated by the symbols on the model lines. The lower 
panel exhibits the [N/q] vs. [a/H] data for DLAs (open circles), 
and upper panel tiie data for dwarf galaxies from van Zee et al. 
(1997a,b) (crosses) and Izotov & Tlman (1999) (plus signs). The 
a-element abundance is given by O for the model predictions and 
for the dwarf galaxy data, and for the DLA data, by O or S, or, 
when not available, by dust-corrected Si. The thick solid lines in 
the upper panel are the statistical fits to the DLA data. 



a large number of systems with near-solar [Si/Fe] ratios, un- 
less we assume unrealistic, extremely low values of v (< 0.1 
Gyr-^). 



4.2 One-zone model for disk galaxies 

The disk models predictions are compared to the observed 

values in DLAs (Figures 9 and 10). There is a great disper- 
sion for the [N/a] ratio, but the observed [N/a] ratios are 
confined to the region between r = 4 and 14 kpc, with a 
significant number of systems around r = 8 kpc. The sta- 
tistical trend lines seem to confirm the behaviour predicted 
by the models: the ratio increases with [a/H] rapidly in low 
metallicities and then becomes fiatter. The plateau at higher 
metallicities is due to primary N production by IMS. The 
comparison between the observed [N/a] with the disk model 
predictions indicates star formation time scales in DLAs sim- 
ilar to those found in the solar neighbourhood. The lower 
bound envelope in the [N/a] vs. [a/H] distribution of the 
DLA data is nearly traced by the evolutionary track at r = 4 
kpc. 

With respect to the [Si/Fe] predictions from the disk 
model, the region r « 8 — 18 kpc accounts for most of the sys- 
tems. There are, however, a few systems with higher [Si/Fe] 
ratios implying a faster SFR than that found in outer re- 
gions of disks. The models also present the same increasing 
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Figure 8. [Si/Fe] vs. [Fe/H] observed in DLAs (open circles) 
compared with one-zone models predictions for dwarf galaxies. 
The data points represent the dust-corrected abundances. The 
solid lines represent the grid of models with u = 0.1, 0.3, 1, 3, 19 
Gyr^^ (bottom to top), and the dotted line a model with u = 3.0 
Gyr~^ and M = 10^ Mq. The thick solid lines in the upper panel 
are the statistical fits. 



Figure 10. [Si/Fe] vs. [Fe/H] observed in DLAs (open circles) 
compared with chemical evolution models. The points represent 
the dust-corrected abundances. The solid lines are the predictions 

of the disk model at r = 2, 4, 8, 14, 18 kpc (top to bottom). The 
thick solid lines in the upper panel are the statistical fits. The 
points with error bars have the same meaning as in Figure 9. 




[c/H] 

Figure 9. [N/a] vs. [o/H] observed in DLAs (open circles) com- 
pared with predictions from models of disk galaxies. The solid 
lines are the predictions of the disk model at r = 2, 4, 8, 14, 18 
kpc (right to left). The thick solid lines in the upper panel are 
the statistical fits. The points with error bars refer to DLAs with 
both [N/a] and [Si/Fe] cor determinations, and are identified by 
the right ascension of the background QSO. 



trend of [Si/Fe] with decreasing [Fe/H] of the data. Since 
external regions (r ^ 8 kpc) are favored as sites for DLAs, 
long star formation time scales are typical for these sys- 
tems. A long and continuous star formation allows the iron 
peak group elements to be produced by SN la, giving rise 
to low values of [Si/Fe]. A slow star formation was also in- 
ferred from the analysis of the [N/a] ratios, but with the 
difference that more central regions (r « 4 — 14 kpc) are 
required by the data. Although the [Si/Fc] ratio apparently 
has a smaller dispersion than the [N/a] ratio, a number of 
systems have values of [Si/Fe] too low or too high to be ex- 
plained as arising in one characteristic portion of a disk. In 
addition, there are a few systems with [Si/Fe] lower than 
that expected even in the outermost regions of a disk, and, 
for a number of DLAs, their low [Si/Fe] and high [Fe/H] 
could be achieved in the outer disk but at times later (t ~ 5 
Gyr or later) than those allowed by the age of Universe at 
the DLA rcdshift. 

Figure 11 presents the evolution the metallicity as a 
function of the H I column density predicted by the disk 
model for the radii r = 2, 4, 8, 14, 18 kpc. This model for 
the Galactic disk predicts, at the outermost radius, metal 
column densities ([Si/H]-|-log N(H I)) that are too low at 
all times. The innermost radius that match the metallicity- 
N(H I) constraints is r = 14 kpc, but only at late times, 
t > 3 Gyr, thus limiting the applicability of this model to 
low rcdshift objects only. At the solar radius and r — 4 kpc, 
the predictions match well the observations. The r = 2 kpc 
zone also satisfy the constraints in Figure 11, but only at 
very early times {t ^ 5 x l(f yr). 
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Figure 11. Metallicity vs. H I column density observed in DLAs 
compared with chemical evolution models. The metallicity is 
given by [Si/H]cor- for the DLAs (filled circles). The solid lines 
are the predictions for disk models at r = 2, 4, 8, 14, 18 kpc 
(left to right). The thick line contour encloses the data accord- 
ing to the general constraints on HI column density for DLAs, 
20 < log N{H I) < 22, and the limits on the metal column den- 
sity of our sample, F([Si/H]cor)=[Si/H]cor + logNiHI) < 20.9 
(upper diagonal) and F([Si/H]cor)> 18.9 (lower diagonal). For 
the sake of clarity, the errors are shown only for the two data 
points with the lowest and highest F([Si/H]cor). 



4.3 The chemodynamical model for dwarf 
spheroids 

The predictions for the [N/a] ratio of the chemodynamical 
model for dwarf spheroids are compared to the observed 
ratios in Figures 12 and 13. One characteristic of the model 
is the occurrence of several short starbursts over the galaxy, 
which lead to rapid excursions of the evolutionary tracks in 
the [N/a] vs. [a/H] plane. Most of the systems are young, 
with ages around 0.1 Gyr or less (three systems, with [N/a]~ 
—1.5, have indicated ages of ~ 0.05 Gyr). The DLAs with 
[N/a]> -0.7 and [a/H]< -1.0 could be older systems {t > 
0.1 Gyr), but the crowding of evolutionary tracks prevents 
us from assigning narrow ranges of ages (and of radii) to 
them. 

The lower bound envelope in the [N/q] vs. [a/H] dis- 
tribution of DLAs, which was reproduced by the one-zone 
dwarf galaxy model and by the disk model, representing an 
upper limit on u, is even more apparent in the chemody- 
namical model with u = 3 Gyr~^, and corresponds to the 
radius r ~ 1 kpc. 

The predicted and observed [Si/Fo] ratios in DLAs are 
shown in Figures 14 and 15. Nearly all [Si/Fe] values of DLAs 
are reproduced in the r ~ 1 — 7 kpc region, which is char- 
acterized by a low to intermediate star formation rate. As 
a consequence of the starbursts, at r ~ 3 kpc, the [Si/Fe] 
ratio varies abruptly allowing the evolutionary track of the 
model to cover most of the data points. Multiple starburst 
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Figure 12. [N/a] vs. [a/H] observed in DLAs (open circles) 
compared to the predictions of chemodynamical model for dwarf 
spheroids with f = 3 Gyr"^. The lines represent the model pre- 
dictions at several radii. The thick solid lines in the upper panel 
are the statistical fits. The points with error bars have the same 
meaning as in Figure 9. 



also happen at r ~ 1 kpc, although less dramatically seen 
in the [Si/Fe] vs. [Fc/H] diagram. The evolutionary tracks 
of the model show for [Fo/H] < —2 a trend of increasing 
[Si/Fe] with decreasing metalicity which may be confirmed 
or not by the no Zn sample, with typically lower metallici- 
ties. More reliable abundance measurements in this range of 
metallicities are needed in order to establish whether there 
is discrepancy between the data and the models. Note that 
a potential discrepancy could be used as a test for the SNe 
II yields. 

The occurrence of multiple starbursts and the peculiar 
tracks in both diagrams is explained by the connection be- 
tween star formation and gas flows. In the dwarf galaxy 
model, there is a first infall episode of primordial gas into 
the dark matter halo potential well which ends with a star- 
burst. This first starburst drives an outflow, which is very 
weak and it is reverted into infall by the large outer gas 
reservoir. During the outflow (except in its beginning) and 
the subsequent infall, the gas density is not high enough to 
allow star formation. Since the infall carries little enriched 
gas from the outer parts of the galaxy, there is a decrease in 
the metallicity of the gas as the infall goes on. Eventually the 
second infall builds up a gas density high enough to trigger a 
second starburst. In the inner regions of the galaxy, the cycle 
restarts a number of times giving rise to a series of bursts of 
star formation. In contrast, in the outer (r ^ 4 kpc) regions 
of the galaxy, the cycle happens only once, and the gas suc- 
ceeds in leaving the galaxy. Only this region presents a true 
galactic wind, which is enriched by SNe la. The plateau 
at low, constant [Si/Fe] and [Fe/H] increasing with time. 
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Figure 13. [N/a] vs. [a/H] observed in DLAs (open circles) 
compared to the predictions of clieniodynamical model for dwarf 
spheroids with u = I Gyr^^. The lines represent the model pre- 
dictions at several radii. The meaning of the line styles is the same 
as in Figure 12. The thick solid lines in the upper panel axe the 
statistical fits. The points with error bars have the same meaning 
as in Figure 9. 



Figure 14. [Si/Fe] vs. [Fe/H] observed in DLAs (open circles) 
compared to the predictions of the chemodynamical model for 
dwarf spheroids with i/ = 3 Gyr~^. The lines represent the model 
predictions at several radii. The meaning of the line styles is the 
same as in Figure 12. The thick solid lines in the upper panel 
are the statistical fits. The points with error bars have the same 
meaning as in Figure 9. The polygon encloses the no Zn sample 
defined in Section 2.3. 



which is reached at t « 1 Gyr in the outer regions, marks 
the onset of the wind driven by SNe la. Therefore, dwarf 
spheroids contribute to the enrichment of the intracluster 
medium mainly in iron and not in a elements (see Friaga 
2000). The galactic wind is not a favourable site for DLAs — 
the low densities and high temperatures (T ~ 10^ K) make 
the cooling time too long for the development of thermal 
instabilities into DLAs — which accounts for the absence of 
DLAs in the galactic wind plateau. 

Figure 16 shows the joint evolution of the H I column 
density and of the metallicity predicted by the chemody- 
namical model for several impact parameters. In this paper, 
we define as HI gas that with T < 2 x 10'' K. During the 
evolution of the galaxy, part of the gas is heated mainly by 
SNe to temperatures higher than this limit. The hotter gas is 
not included in the calculation of the column density. At the 
start of the calculations the chemodynamical model exhibits 
N(HI)= 1.23 X lO'^^ cm^^, through an impact parameter of 
1 kpc (i.e. at r = r^,). The initial infall and, after a while, 
a series of starbursts are visible through excursions in the 
log N(HI) vs. [Si/H] diagram. The starbursts have a more 
dramatic effect in the r = 3 kpc region. In this region, the 
starbusts shift along the strip defined by the observational 
limits on the metal column density, toward higher metallic- 
ities and lower values of N(HI). In the outer regions (r = 7 
and r = 12 kpc), after the early infall, the galactic wind 
prevails and N(HI) is reduced dramatically. Note that at 
r = 12 kpc, differently from at more central impact parame- 
ters, there is no initial increase of N(HI) due to the baryonic 



infall, and N(HI) remains below the lower limit for a DLA. 
Rather, the infall accumulates mass only in the inner regions 
of the galaxy. The outskirts of the galaxy (r ^ 10 kpc) would 
host only sub-DL As and Lyman Limit Systems (LLSs) . (sub- 
DLAs are QALs with N(HI)= 10^^ - 10^° cm'^— this range 
is narrower than that of the original definition of sub-DLAs, 
N(HI)= 10^^ - 2 X 10^" cm-2 [Peroux et al. 2002a], because 
we are considering as DLAs the QALs with N(HI)> 10^° 
cm"2; and LLSs are QLAs with N(HI)> 1.6 x lO^'^ cm"^ 
[Tytler 1982]) For a small impact parameter (r = 1 kpc), 
the model satisfies the simultaneous constraints on H(HI) 
and metal column density of DLAs only during its early 
evolution {t = 0.02 — 0.2 Gyr). At larger impact parame- 
ters, those constraints arc satisfied during a long time span 
(from t ~ 0.03 to > 1 Gyr). We did not try to tune the 
model to exactly match the observed values of N(HI). How- 
ever, only a variation of a factor of two in the masses of 
the models accounts for the whole distribution of DLAs in 
the log N(HI)-metallicity plane, as we see from Figure 16, 
where we show the results for a model with u = 3 Gyr~^ 
and Mq = 5 X 10* M0, and another with 1/ = 1 Gyr~^ and 
Mg = 5 x 10* Mq. 
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Figure 15. [Si/Fc] vs. [Fc/H] observed in DLAs (open circles) 
compared to the predictions of the c-liemodynamical model for 
dwarf spheroids with v = I Gyr^'^. The lines represent the model 
predictions at several radii. The meaning of the line styles is the 
same as in Figure 12. The thick solid lines in the upper panel 
are the statistical fits. The points with error bars have the same 
meaning as in Figure 9. The polygon encloses the no Zn sample 
defined in Section 2.3. 



5 DISCUSSION 

5.1 The abundance ratios 

The presence of dust in DLAs is suggested by the undersolar 
[Cr/Zn] and [Fe/Zn] ratios (Pettini et al. 1994, 1999). This 
behaviour is generally attributed to the fact that Fe and 
Cr are refractory elements while Zn is little depleted, as 
seen in the ISM of the Gala^xy. However, dust depletion is 
not so clearly indicated by the a elements: in most of the 
systems the [S/Si] ratio is close to solar (but there are more 
DLAs with suprasolar [S/Si] than with undersolar [S/Si]). 
The values above solar could be explained by taking into 
account the depletion pattern of the local ISM: while Si is 
depleted, S is nearly undepleted. On the other hand, the 
undersolar values, if real, arc difficult to explain in these 
context. They may represent intrinsic values in DLAs, where 
the nucleosynthesis of S could result of contributions distinct 
from those of the local ISM, or could be explained in terms 
of differences in the ionization corrections for S and S. 

In DLAs, the [Si/Fc] ratios, even after dust-correction, 
reveal the same a-clemcnt enhancement observed in metal- 
poor stars of the Galactic disk. However, there are sys- 
tems with nearly solar [Si/Fe] ratios, even at low metal- 
licities, which is similar to the pattern predicted by dwarf 
galaxy models. This is consistent with a multi-population 
scenario for DLAs. However, the dwarf galaxy chemodynam- 
ical model accounts for nearly all the DLA data, while there 
are a number of systems not explained by the disk model. 
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Figure 16. Metallicity vs. H I column density observed in DLAs 
compared to the predictions of the chemodynamical models for 
dwarf spheroids with Mq = 10^ MQand = 3 Gyr~^ (up- 
per left panel) and i/ = 1 Gyr~^ (upper right panel), with 
Mq = 5 X 10® M0 and v = 2, Gyr~i (lower left panel) and 
with Mq = 2 X 10^ Mq and v = I Gyr^^ (lower right panel). 
The lines represent the model predictions for several impact 
parameters. The meaning of the line styles is the same as in 
Figure 12. The metallicity is given by [Si/H]co7- for the DLAs 
(filled circles). The thick line contour encloses the data accord- 
ing to the general constraints on HI column density for DLAs, 
20 < log N(H I) < 22, and the limits on the metal column den- 
sity of our sample, F([Si/H]cor)=[Si/H]cor + log N{Hr) < 20.9 
(upper diagonal) and F([Si/H]cor)> 18.9 (lower diagonal). For 
the sake of clarity, the errors arc shown only for the two data 
points with the lowest and highest F([Si/H]cor). 



The pure chemical evolution models, both for disk sys- 
tems and dwarf galaxies, present an apparent discrepancy 
(which is more serious for the dwarf galaxy models) be- 
tween the fits to the [Si/Fe] and [N/Q]i ratios. In the disk 
model, the [Si/Fe] ratios are better reproduced at larger radii 
than those required by the [N/a] ratios, and, in the dwarf 
galaxy model, by lower values of ly than those implied by the 
[N/a] ratios. This difference could be related to an under- 
estimated dust correction for Si, but it does not seem to be 
the case, since the trends [Si/Fe]cor- vs. [Fe/H]cor and [S/Zn] 
vs. [Zn/H] are similar. 

Another effect that could contribute to the above dis- 
crepancy is that different samples are used in the com- 
parisons to the models. The DLA sample with [N/a] val- 
ues is smaller than the sample with [Si/Fc]cor values. If a 
multi-population scenario for DLAs applies, some galaxy 
type missing or a few appearances of some rare galaxy 
type in the small number DLA sample with [N/a] values 
could make difficult the comparison between the two sam- 
ples. In addition, the sample with [N/a] values has higher 
redshifts {z 

median — 3.02) than the sample with [Si/Fe]cor 
values (^:median = 2.07), then one cannot exclude the possi- 
bility that the samples are drawn from different populations 
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(different galaxy types or galaxies at different evolutionary 
stages). 

Differently from one-zone models, the chemodynami- 
cal model accounts for nearly all the DLA data. One single 
model (the ly ^ 3 Gyr"\ M = 10^ model) is able to 
reproduce [a/Fe] and [N/q] ratios and HI column densities 
of five out of six DLAs, if we consider a range of radii for the 
DLA sites — from central (r ~ 1 kpc) to relatively outly- 
ing regions (r ~ 7 kpc). Therefore, from the chemical point 
of view, the DLAs could represent a single type of system 
(the dwarf spheroids in our case). On the other hand, direct 
imaging (Le Brun et al. 1997; Pettini et al. 2000) and kine- 
matical data (Prochaska ct al. 1998) indicate that at least 
some DLAs are associated to disks. It is this sort of non- 
chemical evidence which makes a case for a multi-population 
scenario for DLAs. 

We should note that for the system of our sample with 
the highest [Si/Fe], the Zabs = 3.025 DLA toward Q0347-38, 
Lcvshakov et al. (2002) report [Si/Fe]= 0.95 ± 0.02, which 
is significantly higher than the value [Si/Fe]= 1.17 ± 0.03 
of Prochaska et. (2002b) listed in Table 1. The reason for 
this difference is the lower A'^(Si''') measured by Prochaska 
ct. (2002b), and not a higher II(III), because Prochaska et. 
(2002b) adopt log-ff(HI) = 20.626 ± 0.005 from Levshakov 
et al. (2002). Prochaska et. (2002b) suggest that the dis- 
crepancy is due to significantly telluric line-blend in the 
spectrum obtained with the UVES-VLT by Prochaska et. 
(2002b). For this reason, we have adopted the [Si/H] do- 
termination of Levshakov et al. (2002). However, if the high 
[Si/H] is confirmed by further observations, it will pose prob- 
lems both for the dwarf galaxy chemodynamical model and 
for the disk model, because the resulting [Si/Fe]=0.74 will 
be much higher than the predictions of any model presented 
here. If this is the case, it may represent a system with a 
SFR higher (e.g. v = 10 Gyr~^) than that assumed in the 
dwarf galaxy chemodynamical model {u = 3 Gyr~^). 

5.2 The role of gas flows in DLAs 

The failure of the one-zone model for dwarf galaxies in re- 
producing the abundance ratio data is related to a general 
inability of the clEissic wind models, which lead Matteucci et 
al. (1997) to propose a multi-burst model for DLAs, based 
on a single-zone chemical evolution model for dwarf irregu- 
lar, with a series of parameters describing the burst - burst 
duration, number of bursts, interval between bursts. In the 
chemodynamical model, multiple star bursts appear in a 
natural way, and no ad hoc parameters describing them are 
needed. The inner (r ^ 4 kpc) region exhibits a series of 
short bursts of star formation, which can been seen very 
clearly in the [Si/Fe] vs. [Fe/H] plot, and which cover all the 
values derived for the DLAs. The connection between the 
gas fiow evolution and the star formation rate is the reason 
behind the ability of the chemodynamical model in repro- 
ducing the wide range of abundance ratios seen in DLAs. 

The success of the chemodynamical model for dwarf 
galaxies does not represent an argument against the disk 
scenario for DLAs, because we chose a simple model for the 
disk evolution, that with a single infall of primordial gas, 
in order to discern the effect of gas flows in the evolution 
of the system (for the same reason, we studied the classic 
galactic wind model for dwarf galaxies) . The disk paradigm 



for DLAs should be tested by more sophisticated models, 
such as the double infall model for the Galaxy (Chiappini 
et al. 1997). Moreover, the fact that this model considers a 
more elaborated form for infall stress the importance of gas 
flows in chemical evolution of galaxies. Kinomatical data, 
revealing the velocity field underlying the DLAs would be 
very useful for checking distinct DLA scenarios. However, 
as a preliminary conclusion, our results seem to exclude an 
outflow scenario for DLAs, since our model resolves outflows, 
and no DLA is found in the outflow stage/region. 

5.3 The epoch of galcixy formation 

An important implication of the study of the chemical evolu- 
tion of DLAs is to investigate the cosmic progress of galaxy 
formation. The models provide chemical clocks to infer the 
formation epoch of spheroids using the DLA chemical data 
as constraints. We will use as simultaneous constraints on 
the epoch of galaxy formation zap, the [N/q] vs. [a/H] 
and [Si/Fe] vs. [Fe/H] data of six DLAs, which are com- 
pared with the predictions of the chemodynamical models 
with Mg = Mo and ^ = 3 and 1 Gyr"^ The sys- 
tems can be divided in 2 groups according to the region of 
the galaxy reproducing the [N/a] ratio: the inner DLAs, lo- 
cated at r ~ 1 kpc (DLAs toward Q1104-1805, Q2243-6031 
and Q0347-38), and the outer DLAs, in the r ~ 3 — 7 kpc re- 
gion (DLAs toward QOOOO-2620, QOlOO-l-13 and Q1331-M7). 
The model with v = 3 Gyr~^ best fit the inner DLA data. 
The system toward QOOOO-2620 is fitted only by the u = 1 
Gyr"^ model, but it probably has Mq = 2 - 3 x 10^ Mq, 
from the constraints on N(HI). We can see that the inner 
DLAs would be within very young galaxies, with ages ~ 0.1 
Gyr. Only far from the centre, we find older systems: the 
DLAs toward QOOOO-2620, QOlOO+13 and Q1331-1-17, with 
ages of 0.1-0.5, 0.1-0.3 and 0.5-1.0 Gyr, respectively, cor- 
responding to ZGF = 3.558-4.411, 2.392-2.574, and 2.075- 
2.465. The galaxy formation epoch of a DLA in very young 
environment is, of course, only slightly higher than its Zabs- 
Assuming an ago of 0.1 Gyr, the DLAs toward Q0347-38, 
Q1104-1805, and Q2243-6031 have been formed at zgf = 
3.160, 1.711, and 2.414, respectively. These numbers charac- 
terize a typical epoch of formation for these systems around 
ZGF ~ 2 — 3.5 within a range zgf = 1.7 — 4.5. Higher values 
of zgf could be achieved if we consider the systems with 
highest values of [Fe/H], supposed to be the most chemi- 
cally evolved, and good candidates to be old systems. This 
is the case of the DLAs toward Q0551-366 (zats ~ 1.962) 
and Q0201-F36 (Zabs = 2.463), with [Fe/H]cor = -0.013 and 
-0.191 respectively, which are in a region of the [Fe/H] vs. 
[Si/Fe] plane corresponding to galaxy ages of 1-2 Gyr. Very 
early galaxy formation could occur: zgf = 2.791-4.579 and 
3.794-8.147, for Q0551-366 and Q0201+36, respectively. One 
general conclusion is that DLAs constitute a relatively young 
population with ages not exceeding 2 Gyr, or even 1 Gyr. 

5.4 Neutral gas content of DLAs 

While the directly observable baryonic content of galaxies 
at the present epoch is concentrated in stars, in the past, 
this must be in the form of gas. The near coincidence of 
the inferred baryonic mass density of DLAs (Odl^) at ^ ~ 
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Figure 17. Evolution of the baryonic content of tiie models with 
Ma = 10^ Mq, for 1/ = 3 Gyr~i (upper panel) and v =1 Gyr~i 
(lower panel). The thin solid lines give the evolution of the mass 
in neutral gas for three temperature limits (T < lO**, 2 X 10** and 
3 X 10'' K) . The solid line is the evolution of the total amount of gas 
(at any temperature). The dashed line represents the evolution 
of the stellar mass. Also shown the evolution of the neutral gas 
fraction with respect to the stellar mass at z = 0, assuming that 
the galaxy has been formed at ^ = 4.5 (dotted line). 



2 and the mass density in stars today (Storrie-Lombardi 
& Wolfe 2000) makes the DLAs the prime candidates to 
the progenitors of modern galaxies. Since the DLAs arc the 
dominant reservoir of neutral (HI) gas at z > available 
for galaxy formation, it is of interest to establish whether 
the dwarf galaxy chemodynamical model predicts enough HI 
gas to make the DLAs important contributors to the present 
day stellar mass of the universe. In addition, the evolution 
of the H I mass in our models should be compared to recent 
estimates of the evolution of 0,dla (Rao & Turnshek 2000; 
Storrie-Lombardi & Wolfe 2000; Peroux et al. 2001). 

Figure 17 shows the evolution of the HI gas for the 
models with Ma = 10^ M©, and u = 3 Gyr~^ and v = 1 
Gyr~^. We assume that the neutral gas contains only HI and 
no H2, since our calculations do not include molecules. Our 
adopted cooling function is valid down to T = 1900 K, where 
the H2 fraction is still unimportant. Arry gas that cools below 
this temperature is removed fronr the fluid according to our 
prescriptions for thermal instabilities (Friaga 1993, Friaga & 
Jafelice 1999, Gongalves & Friaga 1999), and is used in star 
formation. The identification of the HI gas to the neutral 
gas is supported by the fact that in DLAs A'^(H2) < A''(H I) 
(Petitjcarr, Srianand & Ledoux 2002). We consider as HI 
gas that with T < 2 x 10'' K. This upper temperature limit 
is somewhat higher than the typical ~ 10* K of the warm 
rncdiunr in the Galactic disk, but in the halo warm mcdiunr, 
the temperatures may be considerably higher than in the 
disk warm medium (Sembach & Savage 1996). In addition, 
Prochaska et (2002a) derived from the Doppler parameters 



Figure 18. Comparison of the observed evolution of the Qg due 
to the DLAs to the predictions of the dwarf galaxy chemodynam- 
ical models with Mq = 10^ , assuming that the dwarf galaxy 
population has been formed at z = 4.5. The upper panel shows 
the predictions of the models u = 3 Gyr^' (lower curve) and 
u = I Gyr~' (upper curve), and the lower papel shows the con- 
tribution of a 1/2-1/2 mixture of models with 1/ = 1 and 1/ = 3 
Gyr~^. The open circles at low redshift are the measurements 
from Rao & Turnshek (2000), while the filled circles axe the high 
redshift data from Peroux et al. (2001). Vertical error bars corre- 
spond to l-cr uncertainties and the horizontal error bars indicate 
bin sizes. 

of the lines of the Zabs = 2.625 DLA toward GB1759+7539 
the upper limit T < 3x 10** K. In order to assess the choice of 
an upper temperature limit for the HI gas. Figure 17 shows 
the evolution of the mass in HI gas, assuming three upper 
temperature hmits (T < 10*, 2 x 10* and 3 x 10* K), and also 
the evolution of the total amount of gas, at any temperature. 
All the models are run until t = 13 Gyr. We see that, at 
late times t > 1 Gyr, the high temperature gas makes an 
important contribution to the total gas content of the galaxy. 
Most of the hot gas is confined to the outer regions of the 
galaxy (r ^ 7 kpc) and can achieve temperatures higher 
than 10^ K. In a complete census of metal and mass in gas, 
this hot gas should be included in addition to the HI gas. 
Figure 17 also shows the evolution of the stellar mass, and 
of the HI gas fraction with respect to the stellar mass at 
z — 0, assuming a galaxy formation epoch zgf = 4.5. The 
HI gas fraction allows us to relate the HI gas in DLAs to the 
present day galaxy stellar content. As we see, in the model 
with ly — 1 Gyr~' the baryonic content is dominated by the 
total gas at all times, and the stellar mass becomes larger 
than the HI gas mass only for t ^ 10 Gyr. By contrast, in 
the model with z/ = 3 Gyr~', the gas is quickly consumed 
by star formation, and the fraction of HI gas is low during 
most of the evolution of the model. 

Figure 18 shows the evolution of the cosmic HI gas den- 
sity (in units of the critical density) predicted the models 
with Mg = 10^ Mq, and v = 3 Gyr"^ and 1^ = 1 Gyr"\ 
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compared to the observed floLA evolution. Wc combine the 
HI gas fraction calculated from our models with the recent 
estimate of n,n in stars in the local Universe by Cole et 
al. (2001), n^z = 0) = (0.0037 ± 0.00056) ft^^- We have 
assumed that the models arc typical of the dwarf galaxy 
population, and that population was formed at zqf = 4.5 
and has evolved until z = 0, giving rise to the dwarf galaxy 
contribution to Qt{z = 0), i.e. ildwarf = fl*{M< 0.2M*) = 
0.00098/Ito . where M* = (1.44 ± 0.37) x lO")/!^ is 
the characteristic mass of the z = stellar mass function 
(Cole et al. 2001). The v = I Gyr~^ model is potentially a 
more important contributor to Q,dla, due to its larger HI 
fraction. However, even in the scenario in which each model 
(with u = 3 Gyr^^ and ly = 1 Gyr^'^) contributes to a half 
of ildwarf, the dwarf galaxy population can contribute to a 
half of Qdla aX z > 2, and accounts for nearly all ^dla at 
z > 3.5, although at low redshifts its contribution is less im- 
portant. Therefore, the dwarf galaxies can be an important 
contributor to Qdla, specially at high redshifts. 

It is interesting to note that Boissier, Peroux & Pet- 
tini (2003), using a family of models for disk galaxies, have 
found that also the disk galaxies could account for about 
half of Qdla, however, the contribution of the disk popu- 
lation becomes less significant at high redshift, which is the 
opposite trend of that we find for the dwarf galaxies. This 
complementarity between the two scenarios, again, argue 
for a multi-population scenario for DLAs, in which dwarf 
galaxies arc more important at higher redshifts and disks 
dominate at low redshifts. 
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Figure 19. Mctallicity vs. N(HI) data for several classes of high 
z objects: LAF clouds, LLSs and DLAs, in order of increasing 
N(HI). For DLAs, the metallicity is given by [Si/H]cor derived in 
this paper, and, for LLS and LAF, it is given by [O/H], estimated 
from the observed [C/H] (LAF data taken from Lu 1991 and Lu 
et al. 1998a, and LLS data from Steidel 1990) assuming [0/C] = 
0.3 (Viegas et al. 1999). The curves in the right-upper corner 
are the predictions for LBGs from chemodynamical models with 
several masses, and the metallicity is [O/H], averaged over a 10 
kpc projected radius; the evolution is shown from 0.1 Gyr (lowest 
point of the curve) until 1 Gyr. The thick solid lines are the 
statistical fits to the data. 



5.5 The high redshift cosmic metal content and 
the relation between DLAs and LBGs 

One possible consequence of dust depletion are selection ef- 
fects in optically selected samples. Boisse et al. (1998) iden- 
tified a trend of decreasing mctallicity with increasing N(HI) 
for DLAs. They interpreted the lack of systems with large 
N(HI) and high metallicities as a selection effect due to dust 
obscuration: systems with large N(HI) and high metallic- 
ity (and consequently high dust content) would obscure the 
background QSO, and would not be selected in an optical 
survey. This suggestion needs future confirmation, since, in 
a QSO radio-select survey of DLAs (Ellison et al. 2001b), it 
was found that the bias due to obscuration of QSOs by the 
dusty DLAs seems not to be so important. 

The presence of dust depletion and possible selection 
effects rises the question of a possible underestimate of the 
typical metallicity not only of DLAs but also of other high 
z systems. Such could be the case of LBGs. The success of 
the chemodynamical model for spheroids in reproducing the 
properties of LBGs (Friaga & Terlcvicli, 1999) gives us con- 
fidence in the model prediction that the LBG metal content 
is easily underestimated from a straightforward interpreta- 
tion of the observations. The main spectroscopic properties 
of the LBG population at high z resembles those of star 
forming galaxies (Steidel et al. 1996b), with spectra charac- 
terised by a blue ultraviolet continuum with moderate dust 
extinction, strong interstellar absorption, P Cygni CIV and 
NV lines from massive stars, and weak Lyman a emission 
(Lowenthal et al. 1997, Pettini et al. 1998, 2000b). Ahhough 
the LBGs are very important objects in the study of the 
chemical enrichment of the early universe because they are 



the most obvious star forming regions at high z, they are 
not suited for chemical abundance studies. Despite the high 
quality of their spectra, the large equivalent width of the 
saturated interstellar lines and the high velocity dispersions 
in the ISM of these galaxies make abundance determination 
in LBGs difficult. However, combining observational data 
and modelling allows one to infer typical LBG metallicities 
in the range of 0.3 Zq to Zq (Priaga & Terlevich 1999). 

Another question is whether DLAs are related to the 
LBGs. Mo et al. (1999) proposed that LBGs and DLAs 
form distinct populations: the LBGs would be star form- 
ing galaxies with low angular momentum and short collapse 
time scale, while the DLAs would be the initial stages of evo- 
lution of galaxies with high angular momentum and longer 
collapse times. 

In order to investigate these questions, Figure 19 com- 
pares metallicity vs. N(HI) trends obtained in this work for 
DLAs not only to predictions of the chemodynamical model 
for LBGs (Friaga & Terlevich 1999), but also to data of Lya 
forest (LAF) clouds and LLSs. As we sec from Figure 19, 
the most metal-rich DLAs could be Mq ^ lO" M© LBGs 
in their early evolutionary phases (the 10^^ M© model is the 
progenitor of a present-day ~ 0.3 L* spheroid), but most of 
DLAs are not explained as LBGs. The LBGs typically arc 
sub-L* systems: the characteristic velocity dispersion of the 
emission lines, aem « 70 km"'^ (Pettini et al. 1998; Pettini et 
al. 2001), is reproduced by Mq « 10^° M© models (Friaga 
& Terlevich 1999). However, it seems that most of the DLA 
galaxies would be still fainter than the spectroscopically de- 
tected LBGs. Good fits to the abundance data are obtained 
for the w 10^ M© model. 
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It is interesting to note that other hnes of evidence point 
toward to masses not far from ~ 10^ M© as typical for 
DLAs. A recent Arecibo blind 21 cm survey (Rosenberg & 
Schneider 2003) has revealed that the number of systems per 
unit redshift is dominated by galaxies with HI masses near 
10^ Mq. Zheng & Miralda-Escude (2002) modelled DLAs as 
spherical isothermal gaseous halos ionized by the extralactic 
ionizing background, and concluded that self-shielding im- 
plies a DLA HI mass of « 1.6 x 10* Mq for N(HI)« 3 x 10^° 
cm~^. 

Notice that the normalization of the star formation rate 
for the LEG models is f = W Gyr"'^, whereas for the DLA 
models shown in this paper, v — 3 or 1 Gyr~^. The star for- 
mation proceeds slowly in DLAs, which would give support 
to the scenario for galaxy formation in which there is a di- 
chotomy between the star formation modes for DLAs (long 
star formation time scale) and LBGs (short star formation 
time scale). However, the star formation in DLAs docs not 
need to be much different from that of typical LBGs, since 
a normalization of the SFR only a factor ~ 3 smaller ac- 
counts the properties of a number of DLAs, thus allowing 
for a smoother progression in the evolutionary history of 
DLAs and LBGs rather than a sharp dichotomy between 
the two populations (see M0ller et al. 2002). 

A complete analysis of the metallicity vs. N(HI) relation 
for DLAs should take into account the evolution of metal- 
licity with the redshift, and the fact that this evolution de- 
pends on the column density, and, in particular, that no 
evolution with redshift is observed at log N(HI)> 21 (Boisse 
et al. 1998, Savaglio 2000). On the other hand, for systems 
with low column densities, log N(HI)< 20.8, the metallicity 
does decrease with the redshift (Savaglio 2000). In any case, 
the gaps in Figure 19 highlight domains in column density 
for which further investigation is required. There is some 
suggestion of underestimate of the metal content in DLAs 
from the fact that near log(N(HI))~ 20.5, there seems to 
be a smooth transition between LBGs and DLAs, whereas 
at higher HI column densities, there is a gap between the 
two classes of objects. It is possible that a whole class of ob- 
jects, intermediate between the observed DLAs and LBGs, 
are missing at log(N(HI))<; 21. These could be high metal- 
licity dusty DLA-like objects or LBG-like objects that are 
too faint or have too red G — TZ colours to be classified 
as LBGs. Also on the low N(HI) side of the DLA distri- 
bution it is possible to hide significant amounts of metal. 
The trend of increasing metallicity with decreasing N(HI) 
extends down to the lowest values of N(HI) accepted for 
DLAs (log(N(HI))« 20), and then there is a new gap, filled 
with a few LLSs and sub-DLAs. This gap suggests that a 
large number of these objects is missing. In addition, the ion- 
ization corrections for these systems are very complex, and 
their chemical abundances easily could be underestimated. 
Proper consideration of the radiation fields and adequate 
application of photoionization codes could lead to metallic- 
ity estimates for some LLSs in the range of 0.1 Zq to Zq 
(Viegas & Friaga 1995). The ionization corrections should 
include high ions, which can make a sizeable contribution 
to the total metal column density, specially if we take into 
account that, as said in Section 4.3, these systems are pref- 
erentially located in the outer regions of the gala^xy, where 
the hot gas concentrates and gives rise to a local EUV/soft 
X-ray ionizing radiation field (Viegas & Priaga 1995, Viegas 



et al. 1999). In this sense, the data points of LLSs in Fig- 
ure 19 may represent lower limits to the actual metallicities 
of these systems. Therefore, these low column density sys- 
tems might represent important reservoirs of metals at high 
redshifts (see Pcroux et al. 2001). If we are missing impor- 
tant populations of chemically evolved objects, the overall 
metal content in the Universe at high redshift has been un- 
derestimated, in this way alleviating somewhat the problem 
that the density of metal detected at high redshift is smaller 
than the integrated total metal production associated with 
the star formation activity at high-redshift (Pettini 2000, 
Pagel 2000). 



6 SUMMARY AND CONCLUSIONS 

In this paper we have explored several scenarios for the 
nature of the DLAs by means of a study of the metallic- 
ity evolution of these systems. We have analysed the ob- 
servational data on chemical abundances of DLAs with ro- 
bust statistical methods, and corrected the abundances for 
dust depletion. The results of this analysis have been com- 
pared to predictions of several classes of chemical evolution 
models describing a variety of scenarios for DLAs: one-zone 
dwarf galaxy models, multi-zone disk models, and chemody- 
namical models representing dwarf galaxies. We have found 
that the chemodynamical model for dwarf galaxies is more 
successful than our dwarf galaxy one-zone models and the 
disk models in reproducing the properties — abundances, 
abundance ratios, and HI and metal column densities — of 
DLAs. We have then focused on the dwarf galaxy scenario 
for DLAs, as described by the chemodynamical model, and 
investigated a rmmber of issues: the epoch of formation of 
the host galaxies of DLAs, the contribution of dwarf galax- 
ies toiioLA, the multi-population (disks and dwarf galaxies) 
scenario for DLAs, the relation of DLAs and LBGs, possible 
missing populations of high redshift objects, and "hidden" 
reservoirs of metals at high redshifts. We summarise our re- 
sults as follows. 

1) The a-enhancement observed in DLAs seems to be 
real, even after dust corrections have been applied to the 
data. In our sample, the observed [Si/FeJ^edion = 0.32 is 
lowered to [Si/Fe]median = 0.16 by the dust correction, but 
remains typically suprasolar, a pattern similar to that found 
in metal-poor stars of the Galactic disk. However, there are 
DLAs with nearly solar [Si/Fc] ratio even at low metallici- 
ties, which could be hosted by dwarf gala^xies instead of by 
disks. 

2) The bimodality in the [N/a] vs. [a/H] distribution of 
DLAs pointed out by Prochaska et al. (2002b) is less appar- 
ent in our sample. Instead wc find a lower bound envelope 
in the [N/a] vs. [a/H] distribution, above which the DLAs 
are concentrated, and that corresponds to a local specific 
star formation rate ~ 3 Gyr~^. This envelope is close to 
the evolutionary tracks of the one-zone dwarf galaxy model 
with u = S Gyr^^, and of the disk model at r = 4 kpc, but 
its shape is best reproduced by the dwarf galaxy chemody- 
namical model at r = 1 kpc. 

3) The one-zone models, both for disk systems and 
dwarf galaxies, present an apparent discrepancy between fits 
using [Si/Fe] and [N/a] ratios. In the disk model, the [Si/Fe] 
ratios are better reproduced in larger radii than those re- 
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quired by the [N/a] ratios, and, in the dwarf galaxy model, 
by lower values of f than those implied by the [N/a] ratios. 

4) The chemodynamical model for dwaxf galaxies is 
more successful than the pure chemical evolution models 
in reproducing the properties in DLAs. The connection be- 
tween the gas flow evolution and the star formation rate is 
the reason behind this ability of the chemodynamical model. 
The chemodynamical model presents multiple star bursts, 
fed by gas flows which lead the evolutionary tracks of the 
model to cover all the values exhibited by the [Si/Fe] vs. 
[Fe/H] and [N/a] vs. [a/H] distributions of the DLAs. A 
multi-burst model for DLAs have already be proposed by 
Matteucci et al. (1997), but, in the chemodynamical model, 
the starbursts appear in a natural way, with no need of ad- 
ditional parameters describing them. 

5) Ftom the study of DLAs with simultaneous mear 
surements of [N/a] and [a/Fe] ratios we have inferred the 
ages and formation epoch zgf of DLAs. We have found that 
zgf = 1.7 — 4.5, and ages ranging from 0.1 Gyr to 1 Gyr. 
Higher values of zqf could be achieved if we consider the sys- 
tems with highest values of [Fe/H], supposed to be the most 
chemically evolved, and good candidates to be old systems. 
In this case, we find zgf ^ 3— ^ 8. One general conclusion 
is that DLAs constitute a relatively young population with 
ages not exceeding 2 Gyr, or even 1 Gyr. 

6) The dwarf galaxies cam contribute to about half of 
^DLA &t z > 2, and even more, at z > 3.5, although at 
low redsliifts their contribution is less important. A multi- 
population scenario for DLAs seems to favoured, in which 
the DLA population is dominated by dwarf galaxies at high 
redshifts and by disks at lower redshifts. 

7) The results of the chemodynamical model calls for 
a smoother progression in the evolutionary history of DLAs 
and LBGs rather than a sharp dichotomy between the two 
populations. LBGs and DLAs could constitute a sequence of 
increasing star formation rate, with the LBGs being the sys- 
tems with typically short star formation time scales (--^ 10* 
yr), and the DLAs corresponding to the low v {~ 1 — 3 
Gyr~^) end of the sequence. We also arise the possibility 
that we could be missing a whole population of objects at 
high hydrogen density columns with metallicities intermedi- 
ate between those of DLAs and LBGs. 

8) It is possible that relying only on the observations of 
DLAs could lead to an underestimate of the metal content of 
the high redshift Universe, because we are not properly in- 
cluding in the metal census both the above mentioned high 
column density population and the low column density pop- 
ulation constituted by LLSs and sub-DLAs. 
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